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前 言 


灰色 系统 的 特色 是 研究 “小 样本 ”与 “ 贫 信 息 ” 等 不 确定 性 问题 ， 
因此 通过 充分 开发 .利用 已 占有 的 信息 来 挖 气 系 统 本 身 固有 的 规律 
是 灰色 系统 理论 的 基本 准则 . 我 们 可 以 通过 社会 经济、 生态 等 系统 
的 行为 特征 数据 来 寻求 因素 之 间或 自身 的 变化 规律 . 灰色 系统 理论 
认为 ,尽管 客观 系统 的 表象 复杂 、 数 据 离 乱 ,但 它们 总 有 自身 的 整体 
功能 ,必然 蕴含 着 某 种 规律 ,关键 是 如 何 选择 适当 的 方法 来 挖掘 和 利 
用 . 刘 思 峰 教授 提出 了 冲击 扰动 缓冲 算 子 的 概念 ,并 构造 出 了 一 种 较 
为 广泛 使 用 的 缓冲 算 子 . 一 些 学 者 在 此 基础 上 对 弱化 缓冲 算 子 和 强 
化 缓冲 算 子 进行 了 扩展 研究 . 本 书 在 此 基础 上 ,结合 反 向 累积 法 、 单 
调 函 数 与 新 信息 优先 的 相关 知识 ,构造 了 新 的 缓冲 算 子 ,从 而 推广 了 
缓冲 算 子 的 类 型 ,并 能 更 有 效 地 提高 建 模 预 测 过 程 中 的 精度 , 全 书 主 
要 涉及 以 下 几 个 方面 的 内 容 : 

(1) 在 灰色 系统 理论 缓冲 算 子 公理 体系 下 ,利用 反 向 累积 和 的 概 
念 ,构造 了 一 类 新 的 强 ( 弱 ) 化 缓冲 算 子 ,讨论 了 其 相互 关系 及 其 性 
质 . 并 通过 算 例 验 证 了 该 算 子 序列 的 有 效 性 和 实用 性 ,为 冲击 扰动 系 
统 在 建 模 预 测 过 程 中 出 现 的 定量 预测 结果 与 定性 分 析 结 论 不 符 的 问 


A 


sp 
题 提供 了 解决 的 方法 . 

(2) 利 用 单调 函数 定理 及 新 信息 优先 原则 ,构造 了 新 的 弱 ( 强 ) 化 
缓冲 算 子 ,从 而 大 大 地 拓展 了 弱 ( 强 ) 化 缓冲 算 子 的 应 用 范围 . 对 序列 
前 一 部 分 增长 (衰减 ) 速 度 过 慢 ( 快 ), 而 后 一 部 分 增长 (衰减 ) 速 度 过 
快 ( 慢 ) 的 冲击 扰动 系统 数据 序列 在 建 模 预测 过 程 中 常常 出 现 的 定量 
预测 结果 与 定性 分 析 结 论 不 符 的 问题 ,提供 了 多 种 解决 方案 . 首次 将 
缓冲 算 子 的 构造 与 函数 联系 起 来 ,从 而 为 缓冲 算 子 的 构造 开辟 了 新 

(3) 在 灰色 系统 缓冲 算 子 公理 体系 下 ,证 明了 下 列 结论 : 若 d 是 
由 X(1),…,Xx(n) 所 构成 的 表达 式 , 了 为 严格 单调 递增 ( 减 ) 函 数 ，g 
为 了 的 反 西 数 .在 d 中 ,将 f(zx(k)) 替换 工 (k)(k 二 1,…,n) ,对 得 到 
的 新 表达 式 , 用 函数 g 去 作用 ,最 后 的 表达 式 记 为 e, 若 坟 为 强化 ( 弱 
化 ) 缓 冲 算 子 , 则 e 也 为 强化 (弱化 ) 缓 冲 算 子 . 

(4) 灰 色 问 题 中 的 背景 值 的 构造 和 友 色 GM(1,1) 模 型 病态 问题 
一 直 是 灰色 系统 研究 中 的 两 个 重要 问题 . 利用 连 分 式 与 矩阵 条 件数 
理论 ,对 这 两 个 问题 进行 了 讨论 ,获得 了 一 些 重要 结论 . 

本 书 仓促 之 间 不 免 有 错误 和 不 妥 之 处 ,恳请 专家 、 学 者 和 同仁 多 
加 批评 指正 . 


著 者 
2010 年 5 月 
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1.1 研究 背景 及 研究 意义 


社会 .经济 ,农业 、 工 业 、 生 态 等 许多 系统 ,都 是 根据 研究 对 象 所 属 的 领域 和 
范围 命名 的 ,而 灰色 系统 却 是 按 颜 色 命名 的 . 在 控制 论 中 ,人 们 通常 用 颜色 的 深 
浅 来 形容 信息 的 明确 程度 . 如 艾 什 比 将 内 部 信息 未 知 的 对 象 称 为 “黑箱 ”再 如 
在 政治 生活 中 ,人 民 和 群众 希望 了 解决 策 及 其 形成 过 程 的 有 关 信 息 ,就 提出 要 增 
加 透明 度 . 我 们 用 “ 黑 ” 表 示 信息 未 知 , 用 * 白 ”表示 信息 已 知 ,用 * 灰 ”表示 部 分 信 
息 未 知 ,部 分 信息 已 知 . 相应 的 ,信息 完全 明确 的 系统 称 为 “白色 系统 ”, 信 息 完 
全 未 知 的 系统 称 为 “黑色 系统 ”. 部 分 信息 未 知 、 部 分 信息 已 知 的 系统 称 为 “灰色 
系统 ” 在 灰色 系统 理论 创立 与 发 展 过 程 中 , 邓 聚 龙 教授 发 现 并 提炼 出 灰色 系统 
的 基本 原理 . 读者 不 难看 出 ,这 些 基本 原理 ,具有 十 分 深刻 的 哲学 内 涵 . 

公理 1. 1 (差异 信息 原理 ) 差异 即 信息 , 凡 信息 必 有 差异 . 

公理 1.2 ( 解 的 非 唯一 性 原理 ) 信息 不 完全 ,不 确定 的 解 是 非 唯 一 的 . 

公理 1.3 (最 少 信息 原理 ) 灰色 系统 的 特点 是 充分 开发 ,利用 已 占有 的 最 
少 信息 . 
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公理 1. 4 ( 认 知 根据 原理 ) 信息 是 认 知 的 根据 . 

公理 1.5 (新 信息 优先 原理 ) ”新 信息 对 认 知 的 作用 大 于 老 信息 . 

公理 1.6 ( 灰 性 不 灭 原理 ) 信息 不 完全 是 绝对 的 . 

在 灰色 系统 理论 中 ,由 于 冲击 扰动 系统 的 大 量 存在 ,导致 了 定量 预测 结 
果 与 人 们 直观 的 定性 分 析 结论 大 相 径 庭 的 现象 经 常 发 生 . 问题 的 症结 不 在 于 
模型 的 优 劣 ,而 是 由 于 系统 数据 因 系统 本 身受 到 某 种 冲击 波 的 干扰 而 失真 ， 
因此 寻求 定量 预测 与 定性 分 析 的 耦合 点 是 摆 在 每 一 位 预测 工作 者 面前 的 一 
个 首要 问题 . 

灰色 系统 理论 的 主要 任务 之 一 就 是 根据 社会 ,经 济 、 生 态 等 系统 的 行为 
特征 ,寻求 不 同系 统 变量 之 间或 系统 变量 自身 的 数学 关系 和 变化 规律 ,其 特 
色 是 研究 “小 样本 ”与 “ 贫 信息 ”等 不 确定 性 问题 ,其 中 的 “新 信息 优先 原理 ”是 
灰色 系统 理论 的 信息 观 , 即 认为 新 信息 对 认 知 的 作用 大 于 老 信 息 , 赋 予 新 信 
息 较 大 的 权重 可 以 提高 灰色 建 模 、 灰 色 预 测 、 灰 色 决 策 等 的 功效 ,其 中 的 方法 
体系 一 一 灰色 序列 生成 ,是 指 通过 信息 覆盖 ,选择 适当 的 方法 对 原始 数据 进 
行 挖掘 ,整理 以 寻求 系统 变化 规律 的 技术 . 

刘 思 峰 教授 等 提出 了 冲击 扰动 缓冲 算 子 的 概念 ,构造 出 了 一 种 较为 广泛 
使 用 的 缓冲 算 子 并 研究 了 已 有 算 子 的 关系 及 特性 , 即 当 原始 数据 序列 的 前 半 
部 分 增长 减缓) 速度 较 快 ( 慢 ), 后 半 部 分 增长 (减缓 ) 速 度 较 慢 ( 快 ) 时 , 先 用 
弱 ( 强 ) 化 缓冲 算 子 作用 于 原始 数据 序列 ,然后 利用 GM(1,1) 模 型 进行 预测 ， 
可 以 有 效 地 消除 在 建 模 预测 过 程 中 的 干扰 . 


1.2 基本 概念 


定义 1.1 设 系统 真实 序列 为 
X® = (zr (2)，…zo (za))， 


第 1 章 绪论 


而 观测 到 的 系统 数据 序列 为 
X= (z(1),z(2)，…,zGz)) 
一 (ze@o(1) 十 elyzoO(2) 十 ez (Ca) 十 ev) 
一 X 十 e. 
其 中 e = (el ,ez，…,e,) 称 为 冲击 扰动 项 ,X 称 为 冲击 扰动 序列 . 
一 般 情况 下 ,社会 .经 济 等 冲击 扰动 序列 的 数据 序列 均 为 按时 间 顺 序 记 录 
的 数据 序列 ,最 后 一 项 数据 z(n) 称 为 系统 的 新 信息 . 
定义 1.20 设 系统 数据 序列 为 X 一 (z(1),z(2)，…,z(Gz)) , 称 


xz(b) = 六 xz， 
各 
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Dz = 六 rd) 一 Va, 
=] 4 一 1 tt tml 
DE 立 立 coadon) 

t=l tml 1 = 


= rsp DG + DCt+k— 2)z(), 
i 


分 别 为 1 阶 ,2 阶 ,…,/ 阶 反 向 累积 和 算 子 . 
当 z(1) = x(2) ==… 一 zz) 一 1 时 ,有 


D1= Bant Dnth ol). 


定义 1.3 设 系统 数据 序列 为 X= 二 (zx(1) ,xz(2),…,z(n)), 若 
(1) Vk 二 2,3,… yn ，X(k) 一 zlk 一 1) 之 0 , 则 称 X 为 单调 增长 序列 . 
(2) Vk 二 2,3,…,n，X(k) 一 x(k 一 1) 二 0 , 则 称 X 为 单调 衰减 序列 . 
(3) Ai kz € {2,3,…,n} ,使 得 
zh)— zk m1)>0,rk)— zk —1) 一 0， 
则 称 X 为 振荡 序列 . 


其 中 若 令 M = 加 axz() ，m 一 minz() , 则 称 M 一 m 为 振荡 序列 XX 的 
振幅 . 


图 veaonrmpmma 


定义 1.4 设 X 为 系统 数据 序列 , DD 为 作用 于 X 的 算 子 , X 经 算 子 DD 作用 
后 所 得 到 序列 记 为 XD = (zr(1)d,z(2)d,…,z(n)d) , 则 称 D 为 序列 算 子 . 

对 序列 连续 作用 ,可 得 二 阶 算 子 ,一 直 可 以 作用 到 阶 算 子 , 分 别 记 为 
XD’,…, XD". 

公理 1. 7( 不 动 点 公理 ) 设 X 为 系统 数据 序列 , DD 为 序列 算 子 , 则 有 

znd 一 r(n). 

公理 1.8( 信 息 充 分 利用 公理 ) 系统 数据 序列 X 中 的 每 一 个 数据 工 (k) 
(二 1,2,…,n) ,都 应 充分 地 参与 算 子 作用 的 整个 过 程 . 

公理 1, 9( 解 析 化 与 规范 化 公理 ) 任意 的 x(k)d(k = 二 1,2,…,n) 皆 可 以 由 
一 个 统一 的 z(1) ,zx(2),…,zr(n) 的 初等 表达 式 表达 . 

满足 上 述 三 公理 的 序列 算 子 称 为 缓冲 算 子 , XD 称 为 缓冲 序列 . 

说 明 : 公 理 1.7、1.8、1.9 在 下 文中 简称 为 缓冲 算 子 的 公理 一 、 二 ,三 . 

定义 1.5 设 X 为 系统 数据 序列 , DD 为 序列 算 子 , 当 X 为 单调 增长 序列 . 单 
调 衰减 序列 或 振荡 序列 时 ,者 缓冲 序列 XD 比 数据 序列 X 的 增长 速度 (或 衰减 
速度 ) 增 强 ( 减 弱 ? 或 振幅 减 小 ( 增 大 ), 则 称 缓冲 算 子 DD 为 弱 ( 强 ) 化 算 子 . 

定理 1.1 设 X 为 单调 增长 序列 ,XD 为 缓冲 序列 , 则 

(1)D 为 弱化 缓冲 算 子 全 z() 福 z(E)Q ,k= 二 1,2,… ,ns 

(2)DD 为 强化 缓冲 算 子 3z(k) 宇 z(k)d,k==1,2,…,n 

证 明 : 设 

rh) 一 EEE ,km 1,2,.. 


为 原始 序列 和 中 z(R) 到 zz) ep 


Wad—z(kd , 
rl(k)d = a rr 大 一 1,2，… 


为 缓冲 序列 XD 中 z(A)d 到 z(z)d 的 平均 增长 率 . 由 zx(n)d 一 z(n), 得 


ZX —zxk)_ rmad—r(k)d _ zr(k)d— zr(k) 
A A i pi mn 一 上 十 1 


若 卫 为 弱化 缓冲 算 子 , 则 7r(&) 宇 r(k)d, 即 rk) 一 r(k)d 宇 0, 于 是 z(k) 一 
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ZX(k)d < 0, 即 zA) 去 zk)d, 反之 亦 然 . 

若 也 为 强化 缓冲 算 子 , 则 x(k) 壹 rCk)d, 即 r(k) 一 r(k)d 之 0, 于 是 xz(k) 一 
ZX(R)d 宇 0, 即 zx(A) 三 <(A)d, 反之 亦 然 . 

定理 1.2 设 X 为 单调 衰减 序列 , XD 为 缓冲 序列 , 则 

(1)D 为 弱化 缓冲 算 子 SSz(k) 宇 zr(k)d ,k=1,2,"…,n; 

(2)D 为 强化 缓冲 算 子 ST(k) 坟 TCk)d ,k= 二 1],2,… sn. 

证 明 :与 定理 1. 1 类 似 ,从 上 略 . 

定理 1.3 设 X 为 振荡 序列 ,XD 为 缓冲 序列 . 


(1) 若 DD 为 弱化 绥 冲 算 子 , 则 
max x(k) 之 max{(r(k)d}, min 7(k) < min {xr(k)d}; 
kw lShen 1<ksw TSkt<n 

〈2) 若 万 为 强化 缓冲 算 子 , 则 


max x (k) 过 max{z(k)d}, min x(k) 3 Din {ZCk)d). 
由 定理 1. 1 一 1. 3 可 知 , 单 调 增长 序列 在 弱 ( 强 ) 化 缓冲 算 子 作用 下 ,数据 膨 
胀 ( 葵 缩 ); 单 调 衰减 序列 在 弱 ( 强 ) 化 缓冲 算 子 作用 下 ,数据 萎缩 (膨胀 ). 


1.3 实用 缓冲 算 子 的 构造 


定理 1.4 设 原始 数据 序列 X 二 (zx(1),x(2),…,z(n)), 令 
XD = (x(1)d,z(2)d,*…,r(n)d), 


PD ls dese A 
TI(k)d = 让 下 


则 当 X 为 单调 增长 序列 .单调 衰减 序列 或 振荡 序列 时 ,D 皆 为 弱化 缓冲 算 子 . 
证 明 : 直 接 利用 zx(kt)d 的 定义 , 易 知 定理 成 立 . 
定理 1.5 设 原始 数据 序列 XX 二 (zx(1),z(2),…,z(m)) , 令 
XD = (xr(1)d,z(2)d,*…,r(n)d), 


k= 1,2,0 ,ns 


| 
(Bxemmsmensen 


其 中 


tt 
zh)d = kn 


则 当 XX 为 单调 增长 序列 ,单调 衰减 序列 或 振荡 序列 时 , D 皆 为 弱化 缓冲 算 子 . 
证 明 : 不 妨 设 r(k) 为 单调 增长 序列 , 则 


kr(k)++.…+nr(n) 
zx(k)d—z(k) De ZX(k) 


= K+ Dr(k+1) 一 ze 让 十 … 十 mLzCa) 一 工 )] > 0 
十 "十 n Se 


因此 ， 
Zz(k) < zr(k)d. 
所 以 也 为 弱化 缓冲 算 子 . 
同 理 可 证 , 当 X 为 单调 豪 减 序列 或 振荡 序列 时 ,D 为 弱化 缓冲 算 子 . 
定理 1.6 设 原始 数据 序列 X = (z(1),z(2)，…z(Cz)) , 令 
XD = (x(D)d,r(2)d,,r(n)d), 


其 中 


= Wr(h) twr(n) pe 
zr(k)d 一 直入 ,一 1 2 


则 当 X 为 单调 增长 序列 ,单调 衰减 序列 或 振荡 序列 时 , D 皆 为 弱化 缓冲 算 子 . 
证 明 :不 妨 设 z(k) 为 单调 增长 序列 , 则 


wer(k) 二 十 wor(n) 
Wh 


_ ru)[z(R 十 1) 一 zkR)] 十 … 十 zwn[Lz(z) — x(k)] 之 0 
之 0. 


wh 十 十 to 


ZX(k)d—zr(k)= i 


因此 ， 
Zz(k) < rk)d. 
所 以 D 为 弱化 缓冲 算 子 . 
同 理 可 证 , 当 六 为 单调 衰减 序列 或 振荡 序列 时 , DD 为 弱化 缓冲 算 子 . 
定理 1.7 设 原始 数据 序列 X= 二 (xz(1),z(2),…,X(n)) , 令 
XD = (x(1)d,x(2)d,*…,r(n)d), 


二 | 


其 中 
ZK)d = Lr Ck oe on) J sk = 12 
则 当 XX 为 单调 增长 序列 .单调 衰减 序列 或 振荡 序列 时 , D 皆 为 弱化 缓冲 算 子 . 
证 明 : 不 妨 设 xr(k) 为 单调 增长 序列 , 则 
(Ad 一 [zx%(k) Xe X rm (DT 


三 [zw (RD) X… Xr hk) Js = x(k). 
因此 ， 
zr(k) < rk)d. 
所 以 DD 为 弱化 算 子 . 
同 理 可 证 , 当 X 为 单调 衰减 序列 或 振荡 序列 时 , D 为 弱化 缓冲 算 子 . 
定理 1.8 设 原始 数据 序列 X= (zx(1),z(2),…,zx(n)) , 令 
XD = (zx(Dd,zr(2)d, ,rz(md), 
其 中 


zh 


Zk)d= 3 » I(k),k= 1, ,n, 


则 当 X 为 单调 增长 序列 、 单 调 豪 减 序列 或 振荡 序列 时 , D 箔 为 强化 缓冲 算 子 . 
证 明 : 不 妨 设 z(A) 为 单调 增长 序列 , 则 


(ed 一 zx( 一 开 I » x(k) 一 区 有) 
= LR (zk) — x(n)) So. 
x(n) 


因此 ， 
ZX(k) xk)d. 
所 以 DD 为 强化 缓冲 算 子 . 
同 理 可 证 , 当 X 为 单调 衰减 序列 或 振荡 序列 时 , D 为 强化 缓冲 算 子 . 
定理 1.9 设 原始 数据 序列 X 一 (z(1),zr(2),…,z(n)) , 令 
XD = (z(Dd,z(2)d,…,z(Cz)d)， 


) 灰色 缕 冲 算 子 理论 及 其 


其 中 


(ws 十 … 十 wn)r (kk) 
Wwer(k) 十 十 war (n) 


则 当头 为 单调 增长 序列 、 单 调 衰 减 序 列 或 振荡 序列 时 , D 皆 为 强化 缓冲 算 子 . 
证 明 : 不 妨 设 x(k) 为 单调 增长 序列 , 则 


te et ww) (k) 


TCR)d 一 1,2，…7y 


nT(k) 二- + or Cn) 
Cw) rR) 
rr 
因此 ， 
Zz(k) x(k)d. 
所 以 DD 为 强化 缓冲 算 子 . 


同 理 可 证 , 当 X 为 单调 衰减 序列 或 振荡 序列 时 , D 为 强化 缓冲 算 子 ， 
定理 1.10 设 原 始 数据 序列 六 = (z(1) ,x(2),…,r(n)) , 令 
XD = (x(Dad,r(2)d,e, x(md), 
其 中 


Xx:(k) 
[x (KY oo Cn) ms 


则 当 X 为 单调 增长 序列 .单调 衰减 序列 或 振荡 序列 时 ，D 皆 为 强化 缓冲 算 子 . 
证 明 : 设 x(k) 为 单调 增长 序列 . 则 


zk)d 一 1 2 


Xx(k)d= SR 
[zx (RW or (0 J 
2 
全 [z(em 。 RE = x(k). 
因此 ， 
Zz(k) 之 x(k)d. 
所 以 DD 为 强化 缓冲 算 子 . 


同 理 可 证 , 当 X 为 单调 衰减 序列 或 振荡 序列 时 , D 为 强化 缓冲 算 子 . 


@— 
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1.4 灰色 GM(1,1) 模 型 


定义 1.6 设 
Xeo = (zrO(),r02) ,zn)), 
XY = (zr (2) ,xD Cn)), 
4 
ZV Kk) = > zeo0(i)， 


i=1 


则 称 
ZOR) ar Vk) 一 已 
为 GM(1,1) 模 型 的 原始 形式 . 
定义 1.7 设 
Xe@) = (rH (1) ,rr (2) ,x0 (n)), 
X = (x 1),r D2) ,x (n)), 
4 
TD (RE) 一 zeo(iD; 
i=] 
ZY = (2 (1) ,2 (2) ,2 (n)), 
2 二 二 (zk 一 DD +z Dk)), 
则 称 


TO (RE) Ha dk) 一 
为 GM(1,1) 模 型 的 基本 形式 . 
定理 1.11 设 X'o 为 非 负 序列 ， 
Xe = (xO ,zr 2) ,0 Cn)), 
XY = (xD ,xD 2) ,0 xD (za))， 


‘ 
ZVk) = > zto(i; 
各 


{DD 


) 灵 色 缓冲 算 子 理论 及 其 K 用 


ZY 一 (z0(1)，,z0(2)，…z(z))， 
xD 一 南 (zm 丰 一 D 十 z0(bD) 仆 一 2，vm 
若 5 一 (oa,b)7 为 参数 列 , 且 
Y= (rx®(2), ,7 (zz) )T， 
一 zD(2) 1 
B= 
一 zD(za) 1 


则 GM(1,1) 模 型 zx (Ck) 十 wz (A) = 5 的 最 小 二 乘 估计 参数 列 满足 
和 = (BTB)-BTY. 


定义 1.8 设 
Xe = (x0 x02) ,7 (n)), 
XY = x0 ,D2) ,0 xD Cn)), 
和 
Zk) 一 Dr i); 
j=l 
ZY = (zDD ,2 (2) ,0,2 (n)), 
zO)(A) = 二 (2D +z Ch)), 
则 称 
dz 十 ar 一 
dr to b 
为 GM(1,1) 模 型 
ZV k++ar (k=6 
的 白化 方程 ,也 叫 影子 方程 . 


定义 1.9 设 原始 序列 为 
Xeo = (xO ,x (2) ,0 x0 (zz))， 

相应 的 预测 模型 的 模拟 序列 为 
Xi = (zi (1) ,zx (2) zfo (za) )， 


残 差 序列 为 
eV = (e(l),*,e(n)) 
= (x0 (Dz ,rn) — zx (n)), 


相对 误差 序列 为 
A=( 


则 对 于 和 nn, 称 At 一 


el(l) 
A 


el(k) i x 
油 扣 ;| 为 点 的 模拟 相对 误差, 称 


二 _ 1 
人 A= 和 2 


el(n) 
Ee 


ss 


)= (af 


为 平均 相对 误差. 
定义 1.10 设 Xo = (zo(lD),zo(2)，…yzeo(n)) 为 原始 序列 , 久 0 为 
相应 的 模拟 序列 , e%” 为 残 差 序列 , 则 


元 一 工 > zol,S = 1 Dro) — 2)? 
n t=1 n ta1 
分 别 为 X"” 的 均值 .方差 ; 
E= 1 Dek),S = 1 De(kh) 一 
n t=1 7 4=1 


分 别 为 残 差 的 均值 .方差 ; 
c= 昱 ,2 = p(|e(k) 一 让 < 一 0.6745S)) 


分 别 为 均 方差 比值 与 小 误差 概率 . 


第 名 各 
基于 反 向 累积 法 的 弱化 缓冲 算 子 
理论 研究 


2.1 弱化 缓冲 算 子 序列 


刘 思 峰 ,党 炊 国 等 教授 在 其 专著 (灰色 系统 理论 及 其 应 用 ) 中 构造 了 下 列强 
化 缓冲 算 子 , 设 X 二 《x(1) ,zt(2),…,x(n)) 为 系统 数据 序列 , 令 
XD; = (z(Ddisr(2) di ,zn)di) ,i = 1,2, 
其 中 


ry tt 2 


x (kds = 和 2 


则 当 X 为 单调 增长 序列 .单调 衰减 序列 或 振荡 序列 时 , D， 、D, 皆 为 弱化 组 
冲 算 子 . 

在 此 ,我 们 称 D, 、D: 为 平均 弱化 缓冲 算 子 ,并 利用 反 向 累积 法 ,构建 新 的 
弱化 缓冲 算 子 . 
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定理 2.1 设 X 二 (zx(1),z(2),…,zx(n)) 为 系统 数据 序列 , 令 
XD, = (z(1)d zaod)， 


其 中 


二 二 [Cn 十 1 一 2)z0D0 十 … 十 RE 十 DC 十 1 一 2)z( 有 ] 
CT 
zx(hd, 


it n+ 2) 十 … 十 kk 十 1)…Ck 十 1 一 2)] 


“k=1,2,..,n 
则 当 X 为 单调 增长 序列 ,单调 衰减 序列 或 震荡 序列 时 , D, (! 二 1,2,3,4,…) 为 
! 阶 反 向 累积 和 算术 平均 弱化 缓冲 算 子 (! 一 BCASBO). 其 中 


CD 上 1 一 2) 十 … 十 大 (十 /一 2)] 


二 疡 brn 二 4 一 DD 一 Ch 一 Dk 十 4 一 2)] 


证 明 : 容 易 验证 , 算 子 序列 {D,,l = 1,2,3,…} 满足 缓冲 算 子 的 三 个 公理 ， 
因而 D, 为 缓冲 算 子 . 
(1) 当 X 为 单调 增长 序列 时 , 即 
Z() rk+1) rn), 
所 以 
Xk)di— zr(k) 


Ee 一 ZN 十 … 十 人 十 DR 十 2 一 DCze 十 1) 一 ze)]》 


@ 
i 十 [一 2 十 … 十 有 …(k 十 1 一 2)] 
S70; 
即 
zdi > x(k), 
所 以 DD, 为 弱化 缓冲 算 子 . 
(2) 当 X 为 单调 训 减 序列 时 , 即 


ZX(k) 之 ZE 十 1) 过 … 过 Z(z)， 


灵 色 缓冲 算 子 理论 及 其 几 用 


所 以 
xk)di— zk) 


ae 二 一 2[zGD 一 区 十 … 十 估 十 D… 人 十 2 一 DECz 十 D) 一 ze] 
ri .GOz 十 1 一 2 十 … 十 AR 十 D…(R 十 2 一 2) 
委员 
即 
ZE)di < zk), 
所 以 Di 为 弱化 缓冲 算 子 . 
(3) 当 X 为 振荡 序列 时 , 令 


(CA) = max xi) ,rh) = Pin (i), 


对 Vse {1,2,…,n}), 有 


ZX(k) 之 并 5) rn), 
ZX(h) rs) ,rn), 
则 
ZX(s)di— x(k) 


UD! 3 


{mn 下 全 Dx) z+ tt D's + Dx) 一 区] 


[rtn 十 [一 2) 十 … 十 ss 十 1)*…(s 十 1 一 2)] 


@ si 
<0; 
Zz(s)di— xr(h) 


{z2…(a 十 2 一 2)LzGz) 一 zx(h)] 十 … 十 ss 十 1D*… 


a 
QDI 


—2)[Lz() —zx(h)]} 


让 Do 一 站-… 十 ss 十 1)…(s 十 1 一 2)] 


三 0， 
即 ,对 VsE {1,2,…,n}, 有 
I(Wdi rk) 一 maxz(i), 
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I(sS)di 宇 7(h) = Min Z(i)， 
maxz (Dd < < max (7), 
1 cn 
min I(Ddi 宇 min ID). 
所 以 Di 为 弱化 缓冲 算 子 . 
综 上 可 得 ,{Di,: 一 1,2,3,…} 为 弱化 缓冲 算 子 序 列 . 
当 ! 二 1 时 , 即 为 简单 算术 平均 弱化 缓冲 算 子 


zh) + 十 x(n) 
(Ad ed 


当 ! = 2 时 , 即 为 加 权 算 术 平均 弱化 缓冲 算 子 
好 (CA) 十 … 十 并 CD 


kt +n 

同 理 , 可 以 考虑 三 阶 以 及 更 高 阶 反 向 累积 和 算术 平均 弱化 缓冲 算 子 . 

定理 2.2 设 X= (zx(1),x(2),…,z(n)) 为 非 负 系统 数据 序列 , 令 
XE, = (x(1)e,,***,r(n)e), 


Zz(k)d;, = 


其 中 
(Re = [zm Xe X TRETD JT k= 2, ns 
则 当 X 为 单调 增长 序列 \ 单 调 衰 减 序列 或 振荡 序列 时 , E, (! 二 1,2,3,4,…) 为 / 
阶 反 向 累积 和 几何 平均 弱化 缓冲 算 子 (1 一 BCGASBO). 

证 明 :容易 验证 , 算 子 序列 {E,,! = 1,2,3,…) 满足 缓冲 算 子 的 三 个 公理 ， 
因而 E, 为 缓冲 算 子 . 

(1) 当 X 为 单调 增长 序列 时 , 即 

0<zh) 二 xz( 十 1) AS), 


ze [ (ee Se (0 0 


宇 1， 


x(k)e > rk), 
所 以 EE, 为 弱化 缓冲 算 子 . 
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(2) 当 X 为 单调 衰减 序列 时 , 即 
ZX(k) 宇 x(k 十 DD) 宇 … 宇 Tn) 之 0， 


Ce 区 2) St (2 ) EE 2, 


x(k)e, < zr(k), 
所 以 已 为 弱化 缓冲 算 子 . 
《3) 当 X 为 振荡 序列 时 , 令 
I(k) = maxz(Ci)， 
Ti 
ZX(h) = min zi， 
1Si<n 
对 VsE {1,2,…,n}, 有 
I(k) S75) ,Tn), 
Th) Sr(s) ,rn), 


则 
; wi 9) Fs mm 
人 
和 -的 站 人 To 


这 
即 , 对 VsE {1,2,…,n}, 有 
Z(s)el zr(k)= max zi), 
xX(s)e rh) = min7z(i), 
IE 
Max (Dd < maxzli), 
min x (Dd > min x (i). 
1<ien 1<isn 
所 以 Ei 为 弱化 缓冲 算 子 . 
综 上 可 得 , {E,,l 二 1,2,3,…} 为 弱化 缓冲 算 子 序列 . 


8—— 
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当 = 1 时 , 即 为 简单 几何 平均 弱化 缓冲 算 子 
zk)e: = [zx(k) XX… X rm ; 
当 1 二 2 时 , 即 为 加 权 几 何平 均 弱 化 缓冲 算 子 ， 
zk)es = [x(k) Xe Xr" (JT. 
同 理 , 可 考虑 三 阶 以 及 更 高 阶 反 向 累积 和 几何 平均 弱化 算 子 . 
定理 2.3 设 XX= (x(1),z(2),…,r(n)) 为 非 负 系统 数据 序列 ， 


(有 )ei = [rmt® XX XK TR) DD J 
= TD Cn 二 4 一 2)zx(n) 十 …* 十 k*…(k 十 1 一 2)zx(k)] 
x(k)di ， 
To 十 一 2 十 … 十 和 (十 (一 2)] 
则 


ze SzrR di = 1,2,°°, 
当 且 仅 当 z(1) = x(2) 二 … 二 zx(n), 即 序 列 为 常数 序列 时 ,等 式 才 成 立 . 
证 明 : 由 定理 2. 1 与 定理 2. 2 知 定理 2. 3 成 立 . 


2.2 弱化 缓冲 算 子 序列 的 性 质 


缓冲 算 子 满足 的 不 动 点 公理 体现 了 冲击 扰动 系统 对 其 真实 行为 数据 序列 
新 信息 的 优先 原则 ,对 于 !(! = 1,2,…) 阶 反 向 累积 和 弱化 缓冲 算 子 ,其 对 新 信 
息 zn) 的 作用 还 具有 以 下 规律 : 


性 质 2.1 44 二 1,2,…) 阶 反 向 累积 和 弱化 缓冲 算 子 作用 序列 ,赋予 新 信 


息 zln) 的 权重 


n(n 1—2 


n(n 十 2 一 2) 十 … 十 Rk 十 1)*…(k 十 1 一 2) 
为 有 界 递增 序列 . 
证 明 : 对 YsE {&,k 十 1,…,n}, 有 
开 (& 十 1)…( 开 十 2 一 2) ss 十 1)…(s 十 1 一 2 和 兴 na 十 1)…(a 二 一 2)， 


) 灵 色 缓冲 算 子 理论 及 其 应 用 


故 
(2 十 1 一 2) 十 … 十 和 (十 /一 2) 委 (一 人 十 1)2…(2 十 /一 2)， 


1 im 十 4 一 2) 
Re DEREE TE 


即 权重 数列 为 区 间 [一 二 了 "1 | 上 的 有 界 数列 . 
又 对 !/ 委 /二 1， 


nntl—1) 
TT ED) 


en WNL 2) 


WD ey ,ee 
mn 一 2 二 十 和信 十 全 池 


和 nn/l— 2) 
mntl—2) + tk 2 


即 权 重 数列 关于 /为 单调 增加 数列 . 
性 质 2.2 当 !-> co 时 ， 


nn!— 2) 
nn tl—2) + Tk) 


ws 


证 明 : 由 于 


人 二 De 十 1) + 
kre kt 2) We 


(RHODAktD _ (+D 
RHDRHID +) *™ (7%), 


wn 2 Xn 和 12) 
人 一 1 十 7 一 35 (2 一 1) 


一 co (一 co)， 


Mx 2) 
nn tl—2) + TR TL 2) 


1 
El (一 co)， 


35) 
me(C2 十 /一 区 十 十 1 
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即 当 ! -> se 时 ,灰色 序列 生成 的 第 & 个 缓冲 数据 z(A)d 利用 新 信息 z(z?) 的 权 
重 为 1, 因 此 该 弱化 缓冲 算 子 序 列 充 分 体现 灰色 系统 理论 的 信息 观 . 


2.3 实例 分 析 


【实例 1】 某 市 乡镇 企业 1994 一 1998 年 的 总 产值 (单位 : 亿 ) 分 别 为 
269. 54,284. 45,474. 62,564. 79,597. 6， 
则 每 年 的 增长 率 分 别 为 
31.53% ,23. 68% ,16. 14% ,5. 81%. 

可 见 , 这 个 数据 序列 在 开始 阶段 增长 速度 较 快 ,而 到 后 来 其 增长 速度 变 慢 . 
分 析 原 因 , 可 能 是 开始 阶段 有 相应 的 鼓励 政策 ,使 之 可 以 快速 的 发 展 ,再 就 是 企 
业 发 展 初期 ,由 于 刚刚 起 步 ,基数 较 小 ,开始 进入 的 企业 也 多 是 见效 迅速 的 完全 
竞争 型 企业 ,但 是 随 着 时 间 的 发 展 ,相关 政策 的 收回 和 大 多 数 企业 渐渐 进入 平 
稳 期 ,其 增长 速度 也 就 渐渐 慢 了 下 来 . 

由 此 可 见 ,在 进行 预测 时 ,新 信息 就 代表 了 新 的 情况 和 政策 ,理应 在 预测 时 
有 更 大 的 权重 ,那么 利用 原始 的 数据 序列 直接 建 模 就 是 不 可 取 的 . 为 了 能 够 使 
预测 符合 现实 ,需要 对 该 数据 进行 弱化 处 理 . 下 面 以 1994 一 1997 年 的 数据 为 建 
模 数据 ,1998 年 的 数据 作为 模型 检验 型 数据 . 

为 了 对 比 效果 ,我 们 先 利用 原始 数据 建立 GM(1,1) 的 白化 方程 并 进行 计 
算 , 得 一 步 预 测 结果 如 下 


780. 4701. 
再 利用 Di(! = 1,2,3,…,10) 对 其 作用 ,得 弱化 数据 序列 分 别 为 
XD! = (423. 3500,474. 6200,519. 7050,564. 7900) , 
XD, = (472. 1460.494. 6578,526. 1457,564. 7900) ， 
XD: = (495. 9255,507. 8405,530. 9762,564. 7900) ， 
XD, = (509. 9811,517. 0529,534. 7333,564. 7900) ， 


”一 一 看 
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XDs = (519. 2632,523. 8036,537. 7390,564. 7900)， 
XDe = (525. 8512,528. 9393,540. 1982,564. 7900)， 
XD; = (530. 7701,532. 9653,542. 2475 ,564. 7900) ， 
XDs = (534. 5834,536. 1995,543. 9815,564. 7900)， 
XD, = (537. 6265,538. 8507,545. 4679,564. 7900) ， 
XDio = (540. 1117,541. 0611,546. 7560,564. 7900). 
依次 对 上 述 弱化 序列 建立 GM(1,1) 的 白化 方程 并 进行 计算 ,得 一 步 预测 
结果 分 别 为 
616. 1228 
602. 5447 
594. 1035 
588. 4716 
584. 5028 
581. 5836 
579. 3619 
577. 6233 
576. 2312 
575. 0948 
同样 ,利用 已 (! = 1,2,3) 对 其 作用 ,得 弱化 数据 序列 分 别 为 
XE, = (408. 2480,468. 8396,517. 7457,564. 7900) ， 
XE, = (460. 9872,489. 3111,524. 2184,564. 7900) ， 
XE = (487. 6458,503. 1021,529. 1260,564. 7900). 
依次 对 上 述 弱化 序列 建立 GM(1,1) 的 白化 方程 并 进行 计算 ,得 一 步 预 测 
结果 分 别 为 
620. 1246 
606. 0519 
597. 0748 
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将 如 上 所 得 数据 绘制 成 表 2. 1: 


表 2.1 
弱化 算 子 作用 | 一 步 预测 结果 实际 结果 一 步 预测 误差 

XD， 616. 1228 597. 6000 3.0995% 
XD; 602. 5447 597. 6000 0. 8274% 
XD; 594. 1035 597. 6000 0. 5851% 
XD， 588. 4716 597. 6000 1.5275% 
XD | 584.5028 上 E 597. 6000 2. 1916% 
XDs 581. 5836 597. 6000 2. 6801% 
XD; 579.3619 | 597.6000 3.0519% 
XDs 577. 6233 597. 6000 3. 3428% 
XD。 576. 2312 597. 6000 3.5758% 
XD 575. 0948 597. 6000 3.7659% 
XE 620. 1246 597. 6000 3.7692% 
XE; 606. 0519 597. 6000 1.4143% 

597. 0748 597. 6000 0.0879% 


由 上 表 的 结果 可 得 : 
(1) 算 子 D, 在 1 一 1,2,3 时 一 步 预测 误差 逐渐 减 小 ,在 ! = 3 时 的 一 步 预测 
误差 最 小 ,为 0. 5851%, 从 / 宇 4 直 到 10 又 逐渐 增 大 . 
(2) 算 子 已 在 ! 二 1,2,3 时 一 步 预 测 误差 逐渐 减 小 ,在 /二 3 时 的 一 步 预测 
误差 最 小 ,为 0.0879% ,效果 远 高 于 Ds; . 
(3) 比 起 未 经 弱化 算 子 作用 的 数据 序列 的 预测 值 780. 4701, 其 误差 为 
30. 6008% ,可 见 经 过 弱化 算 子 作用 的 预测 值 要 更 精确 些 . 
【实例 2】 据 统计 ,2002 一 2007 年 中 国 的 高 中 毕业 生 人 数 (单位 :万 ) 分 
别 为 
383. 8,458. 1,546. 9,661. 6,727. 1,788. 3， 
则 每 年 的 增长 率 分 别 为 
19. 36% ,19. 38%,20. 97% ,9. 90% ,8. 42%. 


| | 次 色 缕 冲 算 子 香 认 及 其 应 用 


可 见 前 半 部 分 增长 速度 较 快 ,而 后 半 部 分 增长 速度 较 慢 . 分 析 原因 ,在 2002 
年 左右 的 高 中 毕业 生 , 基 本 上 都 出 生 于 80 年 代 的 人 口 高 峰 时 期 ,而 同时 高 校 也 
在 进行 扩招 ,于 是 带动 了 高 中 学 校 的 扩招 ,因此 在 数据 序列 的 前 几 年 ,毕业 生 的 
数量 有 较 大 的 增长 . 我们 又 知道 ,由 于 国家 的 计划 生育 政策 极 大 地 控制 了 人 口 
的 快速 增长 ,相应 的 ,高 中 毕业 生 的 人 数 开始 渐渐 处 于 小 增长 的 状态 . 所 以 ,如 
果 利 用 原始 数据 直接 建 模 ,数据 信息 的 干扰 性 比较 大 ,而 去 掉 又 有 可 能 减少 相 
关 性 ,为 了 在 新 、 旧 信息 之 间 达 到 一 个 平衡 ,有 必要 用 弱化 缓冲 算 子 对 原始 数据 


序列 作用 . 


为 了 检验 实验 的 效果 , 先 利用 原始 数据 建立 GM(1,1) 的 白化 方程 并 进行 


计算 ,得 一 步 预 测 结果 如 下 
861. 0482. 


再 利用 Di(l = 1,2,3,…,10,…) ,得 弱化 数据 序列 分 别 为 


XD = (555. 5000,598. 4250,645. 2000,694. 
XD, = (614. 8400,631. 3429,660. 2167,697. 
XD, = (644. 6286 ,652. 3000,671. 0935 ,700. 
XD4 = (662. 3671,666. 4043,679. 2231,703. 
XDs = (674. 0540,676. 3760,685. 4708 ,705. 
XDe = (682. 2929,683. 7211,690. 3897,706. 
XD, = (688. 3912,689. 3170,694. 3435 ,708. 
XDs = (693. 0747,693. 7008,697. 5790,709. 
XD, = (696. 7772,697. 2155,700. 2681 ,710. 


XD = (699. 7730,700. 0890,702. 5333 ,711. 
XDa = (713. 5074,713. 5384,714. 0201,717. 
XDu = (720. 3906 ,720. 3931,720. 4704 ,721. 


3500,727. 1000)， 
9889,727. 1000) ， 
9000,727. 1000) ， 
2818,727. 1000)， 
2667,727. 1000) ， 
9462,727. 1000) ， 
3857,727. 1000) ， 
6333,727. 1000) , 
7250,727. 1000) ， 
6882,727. 1000)， 
3963,727. 1000)， 
5255,727. 1000) , 


XDio = (724. 4503,724. 4504,724. 4562,724. 6514,727. 1000) ， 


XD = (724. 6934,724. 6935,724. 6979,724. 8607 ,727. 1000) ， 


22 一 一 一 
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XDix = (724. 9846,724. 9846,724. 9877,725. 1152,727. 1000) ， 
XD = (725. 2130,725. 2130,725. 2152,725. 3177,727. 1000) ， 
XDie = (725. 4505,725. 4505,725. 4519,725. 5311,727. 1000)， 
XDies = (725. 5008,725. 5008,725. 5021,725. 5767,727. 1000) 
依次 对 上 述 弱化 序列 建立 GM(1,1) 的 白化 方程 并 进行 计算 ,得 一 步 预测 
结果 分 别 为 

781. 4491 

764. 1945 

753. 8175 

747. 1908 

742.7321 

739. 5994 

737. 3185 

735. 6079 

734. 2929 

733. 2605 

729. 1229 

727. 6893 

727. 2041 

727. 1866 

727. 1676 

727. 1542 

727. 1417 

727. 1393 

同 理 ,再 利用 Ei(l = 1,2,3) 对 其 作用 ,得 弱化 数据 序列 分 别 为 

XPE, = (540. 7949,589. 2023,640. 7656,693. 5772,727. 1000)， 


三 更 号 
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XE, = (603. 4947,623. 3243,656. 1549,697. 2243,727. 1000)， 
XE, = (636. 1382,645. 6628,667. 4665,700. 1558,727. 1000). 
依次 对 上 述 弱化 序列 建立 GM(1,1) 的 白化 方程 并 进行 计算 ,得 一 步 预 测 


ll 


结果 分 别 为 

786. 5493 

768. 4148 

757. 1667 

将 如 上 所 得 数据 绘制 成 表 2. 2: 
表 2.2 
弱化 算 子 作用 | 一 步 预 测 结果 实际 结果 一 步 预 测 误差 
XD! 781. 4491 788. 3000 0, 8691% 
XD; 764. 1945 788. 3000 3.0579% 
XDs 753.8175 788. 3000 4. 3743% 
XD, 747. 1908 788. 3000 5.2149% 
XD; 742. 7321 788. 3000 5,7805% 
XD， 739. 5994 788. 3000 6.1779% 
XD， 737.3185 | 788.3000 6. 4673% 
XDs 735. 6079 788. 3000 6, 6843% 
XD, 734.2929 | 788.3000 6. 8511% 
XD 733. 2605 788. 3000 6. 9820% 
XD 729. 1229 788. 3000 7. 5069% 
[ XD [ 727. 6893 788. 3000 7. 6888% 
XDio 727. 2041 788. 3000 7.7503% 
| XDin 727. 1866 788. 3000 7.7526% 
XD'1»s 727. 1676 788. 3000 7.7550% 
XDio 727. 1542 788. 3000 7.7567% 
Eo 727. 1417 788. 3000 7.7583% | 

XDies 727. 1393 788. 3000 7.7586% 
XE! 786. 5493 788. 3000 0. 2221% 
XE; 768. 4148 788. 3000 2. 5225% 
757. 1667 788. 3000 3. 9494% 


第 2 章 “基于 反 向 累积 法 的 弱化 缓冲 算 子 理论 研究 这 


从 表 2. 2 可 得 : 

(1) 弱 化 缓冲 算 子 D, 在 ! = 1 时 的 一 步 预测 误差 最 小 ,为 0. 8691%% , 随 着 / 
的 增加 ,其 结果 越 来 越 趋 近 于 z(5) , 即 信息 项 . 

(2) 弱 化 缓冲 算 子 E, 在 ! = 1 时 的 一 步 预测 误差 最 小 ,为 0. 2221%. 

(3) 当 4 = 1,2,3 时 , E, 的 效果 要 好 于 DD, . 

(4) 比 起 未 经 弱化 算 子 作用 的 数据 序列 的 预测 值 861. 0482, 其 误差 为 
9. 2285%, 可 见 经 过 弱化 算 子 作 用 的 预测 值 ,更 加 符合 实际 情形 . 


第 示意 
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3.1 强化 缓冲 算 子 序列 


刘 思 峰 、 党 泡 国 等 教授 在 其 专著 (灰色 系统 理论 及 其 应 用 ) 中 构造 了 下 列强 
化 绥 冲 算 子 , 设 关 = 二 (zx(1) ,zx(2),…,x(n)) 为 系统 行为 数据 序列 , 令 XD 二 (x 
(Ddi,r(2)di,…,x(n)di), 其 中 


一 kt 二 Da 
a yy et 


则 当 X 为 单调 增长 序列 .单调 衰减 序列 或 振荡 序列 时 , D, 皆 为 强化 缓冲 算 子 . 
在 此 ,我 们 称 D, 为 平均 强化 缓冲 算 子 ,并 利用 反 向 累积 法 及 定义 1. 2, 构 建 
新 的 强化 缓冲 算 子 . 
定理 3.1 设 X= (z(1),r(2),…,z(z)) 为 系统 行为 数据 序列 , 令 
XD, = (zx(Ddi,** ,rn di), 


XI(k)di = 


其 中 


lis 十 1)…(z2 十 /一 2) 十 和 (十 /一 2)]z2() 


[To+ 1)…(2 十 1 一 2)z(nz) 十 &…(R 十 2 一 2)z(R) 


zx(k)d 
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则 当 XX 为 单调 增长 序列 ,单调 衰减 序列 或 振荡 序列 时 , D, (! 二 1,2,3,4,…) 为 
! 阶 反 向 累积 和 算术 平均 强化 缓冲 算 子 (! 一 BCASBO) ,其 中 


[nz 十 1)…(m 十 Z 一 2) 十 … 十 AR 十 1)…(R 十 1 2)] 


1 
CT 
一 直 [zOa 十 Da 十 (一 D 一 从 一 DA 十 1 一 2)]， 


证 明 :容易 验证 , 算 子 序列 {D,,/ = 1,2,3,…} 满足 缓冲 算 子 的 三 个 公理 ， 


因而 D, 为 缓冲 算 子 . 
(1) 当 X 为 单调 增长 序列 时 , 即 
Zk) rk+l) Srn), 
所 以 
7 To+ 1)…(m 十 2 一 2) 十 … 十 RE 十 1)…(R 十 2 一 2)]z2(A) 
zx(k)di 一 
CTIC 十 De(n 十 1 一 2)z(n) 十 … 十 Ck 十 1)…(k 十 1 一 2)xz(k)] 
二 Ice+ 1)…(a 十 /一 2) 十 … 十 AR 十 1)…(R 十 2 一 2)]z2(k) 
< Ik). 
re i+ 十 1)…(z 十 1 一 2) 十 … 十 ER 十 1)…( 十 2 一 2)]z() 
即 忆 为 强化 缓冲 算 子 . 
〈2) 当 X 为 单调 衰减 序列 时 , 即 
ZIR) 之 ZI 十 1) 宇 … 之 工 (n)， 
所 以 
区 it “nL 一 2) 十 … 十 RC 十 1 十 1 一 2)Jx2(k) 
x(k)d = 
CI 帮 二 nn 二 DCn 十 1 一 2)z(7) 十 … 十 Ck 十 1)…(k 十 1 一 2)x(k)] 
过 = 十 1)…(2 十 2 一 2) 十 … 十 ACE 十 1)…(E 十 2 一 2)]z2(R) 
(有 ). 
We yi 二 1)…(a 十 2 一 2) 十 … 十 RE 十 1)…(R 十 /一 2)]z(E) 


即 D, 为 强化 缓冲 算 子 . 


和 灰色 缓冲 短 子 理论 及 其 应 用 


(3) 当 X 为 振 功 序 列 时 , 令 


Zz(k) 一 maxZ(i ,x(h) 一 min zi). 


即 
Zz(k) 之 工人 十 1)，…Zz(Cz)， 
I(h) Crh+ ,rn), 
则 
gsi eCn 十 4 一 2) 十 … 十 RC 十 1)*…(k 十 1 一 2)Jz(k) 
Zz(k)di 


一 一 [zz 十 1D)…(n 十 /一 2)z0oD) 十 … 十 RC 十 1)…《R 十 1 一 2)zx(k)] 
U—D! 


万 二 LnCn 十 1)…(n 十 4 一 2) 十 … 十 RCk 十 1)*…(k 十 人 一 2)]x?(k) 
《二 -可 


Ss 二 二 ynt “十 1 一 2) 十 … 十 EC 十 1)…(R 十 1 一 2)]z(k) 
= x(k); 


0 yt "nn 十 4 一 2) 十 … 十 hh 十 1)…(h 十 1 一 2)]z(h) 
(di 一 


(全 vt DD “nL 一 2)z() 十 … 十 hh 十 1)…(h 十 4 一 2)x(h)] 


T= ee *(n 十 1 一 2) 十 … 十 h(h 十 1)*…(h 十 1 一 2)jx2(h) 


a si tb" n 十 1 一 2) 十 … 十 h(h 十 1)*…(h 十 1 一 2)]x(h) 


= zx(h). 
于 是 


TCR)di 三 工 R) = maxz(i), 
1<i 和 nm 


Z(h)di <zr(h) 一 min xi， 
1<i<m 


maxr(i)d, 宇 x(k)d, 宇 rx(k) = max zr), 


1 
min zild <r(h)d,<r(h)= min x). 
所 以 , D, 为 强化 缓冲 算 子 . 


综 上 可 得 , {Di,l = 1,2,3,…} 为 强化 缓冲 算 子 序列 
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当 ! 一 工时 , 即 为 党 氏 强 化 缓冲 算 子 


_ kt Deh) 
A 


同 理 ,可 以 考虑 二 阶 以 及 更 高 阶 反 向 累积 和 算 述 平均 强化 缓冲 算 子 . 
定理 3.2 设 X=(z(1) ,zr(2),…,zx(n)) 为 系统 数据 序列 , 且 zx(i) > 0， 
i= 1 , 令 
XE, = (z(1)e wzGoe)， 
其 中 


Zz(k)e, 一 zh) 
[zr O00 NE x. XK 工 (kN = 


] TREE 


k=1,2,.,n, 
则 当 XX 为 单调 增长 序列 ,单调 衰减 序列 或 振荡 序列 时 , E,(1 = 1,2,3,4,…) 为 / 
阶 反 向 累积 和 几何 平均 强化 缓冲 算 子 (! 一 BCGASBO). 
证 明 :容易 验证 , 算 子 序列 {E,,l = 1,2,3,…} 满足 缓冲 算 子 的 三 个 公理 ， 
因而 已 为 缓冲 算 子 . 
01) 当 X 为 单调 增长 序列 时 , 即 
0<z(k) Crk+1) < rn), 


则 
2 
5 [ET ,。 X 工 Sa Tea 和 nc 

Zz: (k) 
ry Lat J tr 
= 

由 定理 1. 1 一 1. 3 知 , E, 为 强化 缓冲 算 子 . 

(2) 当 X 为 单调 衰减 序列 时 , 即 


Z(R) 之 TI 十 1) 过 … 二 Zoo > 0, 


灰色 组 冲 算 子 理论 及 其 K 用 


则 
CR)ei 本 kD i 
[x (mr 生计 全 XXX 工 (有 人 上 全 ] war 
2 (k) 
之 所 (EPE x Xr kN Jr ns 
= x(k). 
由 定理 1. 1 一 1. 3 知 , E, 为 强化 缓冲 算 子 . 
(3) 当 X 为 振荡 序列 时 , 令 
ZX(k) = maxr (i), 
lSien 
ZX(h) = minz(i)， 
lSics 
即 
Tk) xR), rn) > 0,0<rh) rh+l) ,rn), 
则 
[zx NE x .x x (kN J 
> x*(k) 
[zr (NE XX x (RD ] Am 
= x(k); 
x(h)e= 2 
[x 0) EE? x x x (hE J 
zx (h) 
< 《有 JE 全 XXxr Re Te Te OFDTOTE 
= zx(h). 
则 


ZX(k)e x(k) = maxr(i), 
licn 
xh)e <zr(h) 一 Bin zi), 
maxZ(i)e 之 工人 )e > zx(k) 一 max zt(i), 
1<i<n Sin 
min zx(De <zr(h)e <zr(h) = min 工 (GD). 


1<i<n 


(bd 人 = 5 
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由 定理 1. 1 一 1. 3 知 , E, 为 强化 缓冲 算 子 . 
综 上 可 得 , {E,,! = 1,2,3,…} 为 强化 缓冲 算 子 序列 . 
当 4 = 1 时 , 即 为 党 氏 强化 缓冲 算 子 


Zz(k)e 开外 
[zk) X zx( 十 1) XK" Xr 


同 理 ,可 以 考虑 二 阶 以 及 更 高 阶 反 向 累积 和 几何 平均 强化 缓冲 算 子 . 
定理 3.3 设 X= (z(1),z(2),…,z(z)) 为 非 负 系统 数据 序列 ， 


(De a 
A 人 XX kh) J rn" 


克 一 [n(n 十 Dn 十 1 一 2) 十 … 十 Ck 十 1)*…(kR 十 1 一 2)]z(k) 


it (nn 十 1 一 2)zbn) 十 … 十 kk 十 1)*…* .十 1 一 2)z(D] 


ZC)el < rk)di,l 一 1,2…， 


当 且 仅 当 x(1) == z(2) 一 … 一 z(n) ， 即 为 常数 序列 时 ,等 式 才 成 立 . 


证 明 :由 定理 3. 1 与 定理 3. 2 知 定理 3. 3 成 立 . 


3.2 实例 分 析 


【实例 1】 下 面 以 文献 [1] 中 的 1995 一 2003 年 南京 市 农林 牧 渔 总 产值 数据 


为 例 来 说 明 本 节 提出 的 强化 缓冲 算 子 序列 , 见 表 3. 1. 


表 3.1 1995 一 2003 年 南京 市 农林 牧 渔 总 产值 (单位 : 亿 元 ) 


年 份 1995 1996 1997 1998 1999 

总 产值 86. 1946 91. 9895 94. 2439 96.9644 | 98.9199 | 
年 份 | 2000 2001 2002 2003 

总 产值 | 106.341 | 113.290 | 120.010 | 132.590 | 


, 灰色 缓冲 算 子 理论 及 其 应 用 


以 1995 一 1999 年 数据 作为 系统 行为 数据 ,分 别 在 无 缓冲 作用 、1 阶 .2 阶 反 向 
累积 和 算术 平均 强化 缓冲 算 子 以 及 几何 平均 强化 缓冲 算 子 作用 下 ,建立 GM(1,1) 
模型 , 见 表 3, 2, 模 拟 2000 一 2003 年 农林 牧 渔 总 产值 以 及 误差 , 见 表 3. 2. 

表 3.2 不 同 强化 缓冲 算 子 作用 后 预测 值 及 相对 误差 表 


年 份 2000 年 2001 年 2002 年 2003 年 
无 缓冲 算 子 101. 544 104. 072 106. 663 109. 319 


相对 误差 一 4.51% | 一 8.14% | 一 11.12% | —17.55% 
1 一 BCASBO 102. 94015 | 106. 86788 | 110. 94548 | 115， 17865 | 
相对 误差 一 3.20% | 一 5.67% | 一 7.55% | —13.13% 
2 一 BCASBO 103. 31578 | 107. 56404 | 111.98699 | 116.5918 
相对 误差 | 8 一 5.05% | 一 6.69% | —12.07% 
3 一 BCASBO 103. 54764 | 108. 00397 | 112. 65209 | 117.50024 
相对 误差 —2.63% | ~4.67% | —6.13% | —11.38% 
1 一 BCGASBO | 102.9238 | 106.83738 | 110.89977 115. 11664 | 
相对 误差 一 3.21% | —5.70% | 一 7.59% | 一 13.18% 
2 一 BCGASBO | 103. 30171 | 107. 53742 | 111. 94681 | 116.537 
相对 误差 一 2.86% | 一 5.08% | —6.72% | —12.11% 
3 一 BCGASBO 103. 53611 | 107.9818 | 112.61838 | 117.45404 
相对 误差 —2.64% | —4.69% | 6.16% | —11.42% | 
由 表 3. 2 可 得 : 


(1) 在 没有 进行 缓冲 处 理 情况 下 ,预测 值 的 平均 相对 误差 为 10. 33%6 ,采用 
1 阶 .2 阶 ,3 阶 反 向 累积 和 算术 平均 强化 缓冲 算 子 以 及 几何 平均 强化 缓冲 算 子 
处 理 后 ,预测 值 的 平均 相对 误差 均 比 未 处 理 的 小 . 
(2) 随 着 阶 数 的 增长 ,相对 误差 呈 下 降 的 趋势 . 
【实例 2】 据 统 计 ,1996~1999 年 国民 总 收入 数据 (单位 : 亿 元 ) 序 列 为 
X = (70142. 4916,78060. 8315,83024. 2798,88479. 1548)， 
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则 平均 每 年 的 增长 率 仅 为 8. 1%, 而 1997 年 到 1998 年 的 增长 率 为 6. 4%， 
1998 年 到 1999 年 的 增长 率 为 6. 6%. 这 与 国家 实行 的 五 年 计划 要 求 明显 不 
符合 ,长 期 发 展 下 去 , 必 将 导致 产业 结构 发 展 不 平衡 ,影响 国民 经 济 的 可 持续 
增长 . 因此 ,从 2000 年 开始 ,国家 开始 逐步 调整 产业 结构 ,使 这 种 缓慢 增长 的 
状况 得 到 改善 ,为 了 能 够 及 时 准确 地 把 握 经 济 发 展 趋势 ,对 经 济 的 发 展 作 科 
学 的 预测 ,必须 对 缓慢 增长 的 数据 加 以 处 理 ,使 其 符合 今后 的 发 展 趋势 ,在 此 
基础 上 进行 合理 的 预测 ， 

以 1996 一 1999 年 的 数据 作为 系统 行为 数据 ,分 别 在 无 缓冲 作用 、1 阶 、2 阶 
反 向 累积 和 算术 平均 强化 缓冲 算 子 下 ,建立 GM(1,1) 模 型 ,预测 2000-~2004 年 
的 总 产值 及 其 误差 . 

(1) 无 缓冲 作用 下 : 

原始 序列 为 

X = (70142. 4916,78060. 8315,83024. 2798,88479. 1548). 

GM(1,1) 模 型 为 


i 
内 一 0.6266XO' = 71191. 666486. 
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时 间 响 应 式 为 
Nv (十 1) = 1206295. 830314eaezst — 1136153. 338714. 
由 此 得 ,模拟 序列 为 


A 
X= (70142. 4916,78005. 205721,83049. 417901,88419. 81442). 


残 差 序列 为 
et = (0, 一 55. 629279,25. 138101, 一 59. 34038). 
相对 误差 序列 为 
A 一 (0, 一 0.00071264,0. 00030278, 一 0. 00067067). 
平均 相对 误差 为 


A = 0.056203%. 


) 次 色 缕 冲 牌子 畦 ip 及 其 应 用 


预测 值 为 
X! = (94137. 488013,100224. 895375,106705. 945367, 113606. 092919， 
120952. 439003). 
实际 值 为 
X’ = (98000. 454308,108068. 220558,119095. 689274,135173. 976149， 
159586. 747917). 
预测 值 和 实际 值 的 相对 误差 为 
A' = (一 0.0394, 一 0. 0726, 一 0. 1040, 一 0. 1596, 一 0. 2421). 
(2)1 阶 反 向 累积 和 算术 平均 强化 缓冲 算 子 作用 下 : 
原始 序列 为 
X = (70142. 4916,78060. 8315,83024. 2798,88479. 1548)， 
1 阶 缓冲 序列 为 
XD, 一 (z(l)di ,z(2) 必 ,z(3)di yz(4)cdi) 
一 (61556. 0229,73249, 5903,80383. 5918,88479. 1548). 
GM(1,1) 模 型 为 


dX 一 0.094472Xo) = 63932. 007587. 


dt 
时 间 响应 式 为 
和 (十 1) 一 738287. 535747en ?72 一 676731. 512574. 
由 此 得 ,模拟 序列 为 
入 一 (61556. 0229,73148. 13319 ,80395. 51193,88360. 947255), 
残 差 序列 为 
es 一 〈0, 一 101.45711,11. 92013, 一 118. 207545). 
相对 误差 序列 为 
4 一 (0, 一 0.001385,0. 000148 , 一 0.001336). 
平均 相对 误差 为 


A = 0.095646%. 
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预测 值 为 
X? = (97115. 582851,106737. 611191,117312. 971907， 


128936. 119369,141710. 866306). 
实际 值 为 
X' = (98000. 454308,108068. 220558,119095. 689274， 
135173. 976149,159586. 747917). 
预测 值 和 实际 值 的 相对 误差 为 
A’ = (一 0.0090, 一 0.0123, 一 0.0150, 一 0. 0461, 一 0. 1120). 
(3)2 阶 反 向 累积 和 算术 平均 强化 缓冲 算 子 作用 下 : 
原始 序列 为 
X = (70142. 4916,78060. 8315 ,83024. 2798,88479. 1548). 
2 阶 缓冲 序列 为 
XD;,= (z(1)d:,z(2)d:,z(3)dz ,x(4)d;,) 
= (59330. 0897,72244. 2856 ,80020. 0013,88479. 1548). 


GM(1,1) 模 型 
dx 
一 0.101212X' = 62607. 573071. 
时 间 响 应 式 为 
Ne (p+ 1) = 67909. 266371er or 一 618579. 176671. 
由 此 得 ,模拟 序列 为 
从 =- (59330. 0897,72204. 851743,79895. 468871,88405. 221977)， 
残 差 序列 为 
ee = (0, 一 39. 433857, 一 124. 532429, 一 73. 932823). 
相对 误差 序列 为 
A 一 (0,0. 00054584,0. 00155626 ,0. 00083559). 
平均 相对 误差 为 


A = 0.097923%. 
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预测 值 为 

Xs = (97821. 358122,108240. 417148,119769. 221459,132525. 97123， 
146641. 456266). 

实际 值 为 

X' = (98000. 454308,108068. 220558,119095. 689274,135173. 976149， 
159586. 747917). 

预测 值 和 实际 值 的 相对 误差 为 
A’ = (一 0.0018,0. 0016,0. 0056, 一 0. 0196, 一 0. 0811). 

(4) 将 以 上 预测 值 及 其 相对 误差 做 成 表格 如 下 ; 

表 3.3 不 同 强化 缓冲 算 子 作用 后 预测 值 及 其 相对 误差 
年 份 
2001 2002 2003 2004 


算 子 及 误差 
2000 | 


| 无 缓冲 算 子 | 94137. 488013 [100224. 895375|106705. 945367|113606. 092919| 120952. 439003| 


| 相对 误差 /% —3.94 | 一 7.26 一 10.40 一 15.96 一 24.21 
| 1 一 BCcAspO 97115. 582851 | 106737. 611191|117312. 971907| 128936. 119369|141710. 866306| 
| 相对 误差/% 一 0.90 一 .23 一 1. 50 一 4 61 =i, 00 


|2—BcAsEO 97821. 358122 |108240. 417148 119769. 221459| 132525. 97123 |146641. 456266| 


0.56 一 1.96 


由 表 3. 3 可 得 : 

(1) 在 没有 进行 缓冲 处 理 的 情况 下 ,预测 值 的 平均 相对 误差 为 12. 354%; 采 
用 1 阶 .2 阶 反 向 累积 和 算术 平均 强化 缓冲 算 子 处 理 后 ,预测 值 的 平均 相对 误差 
比 未 处 理 的 小 . 

(2) 随 着 阶 数 的 增加 ,相对 误差 呈 下 降 趋势 . 

(3) 应 用 强化 缓冲 算 子 作 用 后 的 数据 建 模 能 取得 良好 的 预测 效果 . 
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4.1 弱化 缓冲 算 子 序列 


刘 思 峰 ,党 耀 国 等 教授 在 其 专著 (灰色 系统 理论 及 其 应 用 》 中 构造 了 下 列 弱 
化 缓冲 算 子 , 设 X = (x(1) ,zr(2),…,zx(n)) 为 系统 行为 数据 序列 , 令 
XD, = (z(1)di pz(Go)di) ,i = 1,2. 
其 中 
Wer(k) 十 … 十 wrz《n) 


OMT tt 


zk ds = [7 (RW Xe XI MY J ,k= 12 n, 
则 当 久 为 单调 增长 序列 ,单调 衰减 序列 或 振荡 序列 时 , D，、D; 缘 为 弱化 缓冲 算 


和 


在 此 ,我 们 在 弱化 缓冲 算 子 D, 、D; 的 基础 上 ,利用 单调 函数 理论 构建 新 的 
弱化 缓冲 算 子 . 

定理 4.1 设 X==(z(1),zx(2),…,zx(n)) 为 系统 行为 数据 序列 ,是 xz(i)>>0， 
了 >>0,w 二 0. 权重 向 量 zw 一 (zw ,…z)， 太 为 严格 单调 递增 (或 递减 ) 函数 ， 
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全 


8 为 其 反 函 数 , 令 
XDas = (z(1)ds，…zCa)das)， 


其 中 


9 


zh)ds g(t ee 1,2 


则 当 XX 为 单调 增长 序列 .单调 衰减 序列 或 振荡 序列 时 , D; 为 弱化 缓冲 算 子 . 
证 明 : 容 易 验 证 
zds 一 [ee j= xz(n), 


即 D; 满足 缓冲 算 子 公理 一 . 

至 于 缓冲 算 子 公理 二 .公理 三 ,Di 显然 满足 . 因而 , D; 为 缓冲 算 子 . 

假设 7 为 严格 单调 递增 函数 . 

(1) 当 X 为 单调 增长 序列 时 , 即 

0<~<zx(k) rn), 
则 
0~<f(zr(k)) f(rn)), 
0< wt wf TKR) SC wf (zk)) tT Tw f (rn)), 


Ze 十- “+wf(z(n)), 
WW 十 "十 ww 


zds = 8 (AE TE > gf))) = zh). 


Wi 十 十 Ww 


0<f(r(k)) < 


所 以 , D; 为 弱化 缓冲 算 子 . 
(2) 当 X 为 单调 衰减 序列 时 , 即 
T(R) 之 … 之 Z) >0, 
则 
f(r(k)) "f(r(n)) 二 0， 
Owf rk) tw f (Xn)) SC wt ti) Fr)), 


0< we) + 二 < CzC6)， 
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zds — g (LEE tt EY) < g(Cf(z( 昌 )) = zh). 


wh ~ 


所 以 , D; 为 弱化 缓冲 算 子 . 


(3) 当 X 为 振荡 序列 时 , 令 
Z(A) = max zr(i) ,zh) = min Il<i<vZ(i)， 
lSicn 


对 任意 的 i€ {1,2,…,n}, 有 
Zk) Sz ,es TM zh) Sz) ,ez 


f(z(k)) fr) ,f(r(n)) > 0, 
0<fl(rh)) < f(r) , ,flr(n)), 
0<wf z+ twf (rn) < wt tw) fzk)), 
0< (wt tw )f rh)) Swf (zr)) + 二 wf (rn)). 


(ZO7)) 十 … 十 tw CzGD) pz(h)), 


Ti 
0 wi 十 十 ww 


zdy 一 8 2 人 GD 十 全 十 mo 0 fz(k))) 一 (Di 
Wwif (Xz)) + 二 wf (rn)) 


0< f(z(h)) < TT 


zd = (SL ED tf E00 | > go fer))) = zth). 


max (i) > maxz(D)d;, min xz(i) < min z(i)ds. 
lSien leien lei<n lSiSn 
由 定理 1. 1 一 1. 3 知 , D; 为 弱化 缓冲 算 子 . 
同 理 可 证 , 当 f 为 严格 单调 递减 函数 时 , Ds 也 为 弱化 缓冲 算 子 . 
定理 4.2 设 X 二 (x(1),z(2),…,z(n)) 为 系统 数据 序列 , 且 zx(i) > 0， 
了 二 0 二 0. 权 重 向 量 世 一 (zw，…zo)， 太 为 严格 单调 递增 (或 递减 ) 函 数 ， 
8 为 其 反 函 数 . 令 
zd = g{Lf% zk) X…X PCz(D Tr } ,k= 2),n, 


其 中 
XD, = (z(Dds,*%* ,zn)ds), 
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则 当 X 为 单调 增长 序列 ,单调 衰减 序列 或 振荡 序列 时 ，D, 为 弱化 缓冲 算 子 . 
证 明 : 容 易 验 证 
zds = g{[Lf™ (zn) BB}= zm), 
即 D, 满足 缓冲 算 子 公理 一 . 
至 于 缓冲 算 子 公理 二 ,公理 三 ,D, 显然 满足 . 因而 , D, 为 缓冲 算 子 . 
假设 f 为 严格 单调 递增 函数 . 
(1) 当 X 为 单调 增长 序列 时 , 即 
0<zh SZ, 
则 
0 一 ze)) < < frm)), 
0 一 jz(b) 二 [Pa (xh) X…X f(x)) J, 
0<zh) 一 grz(D)) 二 red 一 8{CFa Cz) XX f(x)) ams }, 
所 以 , D, 为 弱化 缓冲 算 子 . 
《2) 当 X 为 单调 衰减 序列 时 , 即 
(和 ) 宇 … 之 zlm) 二 0， 
则 
f(z(k) fr) >0, 
fz(R)) 三 [Fa (zk) Xe Xf (zm)) J > 0, 
Zk) = gf(z(h)) 二 BEC (zk) Xe Xf x) Rs )} = zd >0, 
所 以 , D, 为 弱化 缓冲 算 子 . 
(3) 当 X 为 振荡 序列 时 , 令 
x(k)= Mi) rhy 到 mina), 
对 任意 的 iE {1,2,…,n), 有 
Zk) FID ss Tzh) SI), zn). 
fz CR) 2 fr) ,7s fx) > 0, 
0 < f(z) fz) ,es frm)), 
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0 二 [As (CzGD) XX fr) Ts < f(x)). 


0 一 xz)d = g{Lf% (zDD) XX fs zm rm} gflz(h))) = x(h), 
[Ff (zDD) XX fe zDD) Tm SF f(z)) > 0, 
rd = g{[f® (rz) XX f(z) J 2 gf(zh))) = zx(h) > 0. 
max (i) > maxr(Dd:, minz(i) < min rx(Dd. 
1&i<n 1<i<m lSi<n laicn 
由 定理 1. 1 一 1. 3 知 , D, 为 弱化 缓冲 算 子 . 
同 理 可 证 , 当 f 为 严格 单调 递减 函数 时 , D; 也 为 弱化 缓冲 算 子 
定理 4.3 设 X= (z(1),z(2),…,z(z)) 为 非 负 的 系统 行为 数据 序列 , 且 


Zi) 之 0,f 之 0ww 之 0,wf>>0. 权 重 向 量 分 别 为 w= (ray… yw) 和 w == (wf …， 
由) ,其 中 几 天 0 天 0 着 守之 考 ， 即 


/ 天 
WW WW iT 一 2 


下 为 严格 单调 递增 (或 递 碱 函数, g 为 其 反 函 数 , 则 : 


(1) 当 X 为 单调 增长 序列 时 ,有 
wzGD) Df (x0)) 
5| 二 一 >g| 3 
Dw Pw’ 


(2) 当 XX 为 单调 递减 序列 时 ,有 
Fwif Cz) Dwif ert)) 
SE ， 
Dw } of 
证 明 :假设 f 为 严格 单调 递增 函数 , 则 g 也 为 严格 单调 递增 函数 . 欲 证 


Duwif (xi)) Dw'f (zr)) 
二 | ， 


Sw, wi 
i=k 


i=k 


,对 灰色 缓冲 算 子 理论 及 其 应 用 加 


则 只 需 证 明 
wdcG)) 六 uirczGD) 
站 光一 ， 
Dw Dw 
考虑 Fw Dwif zi)) 一 Pw, Dwif lz) 的 符号 ， 
i=k i=k i=k i=k 
经 计算 ,有 
Fw’ Swf zl)) 一 Dw, Df Cz)) 
i=k j=k i=k i=k 


= [wwaf (zk)) + ww FFCz(R 十 1)) 十 … + wirw,f (rn)) 
wirwef (zk)) + + winwaf (rk)) +wrnwan f(z(k+ 1)) 
+ wwaf (zm) — Lwinwtf (zk)) + ww tn f(rCk+ 1)) 
十 … 十 rotrosf(z(za)) 二 wenwif (zr(k)) 十 …… 十 runrokFCzCR)) 
十 zwrotHiCz( 十 1)) 十 … 十 zewrosFCz(z)) 


> > ww ww) fz) ~ fz))). (x%) 
i=k jetl 
由 于 
, 
Wh i 区 
i 3 
有 
w > wiX YL ,wa 宇 w x 2 
Wi Wiz 
故 
wwXE www, 
Wee Wis wi 
即 


(ww — ww’) 0,7 > i 
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(1) 当 X 为 单调 增长 序列 时 ,有 
(CzG) 一 zGD) 之 0 二 证 
又 因 f 为 严格 单调 递增 函数 , 则 
(f(x) — f(z) 20j >i. 


由 (* ) 式 知 ， 
> py (ww — ww ) (fr) — flr))) 三 0， 
故 
Duwif x0) Pw’f Cz) 
E73 四 > 了 
Dw wi 
jk i=k 
又 因 g 为 严格 单调 递增 函数 , 则 
Dwif lz) Dwif Cz0)) 
br gl ; 
Dw Dw 
(2) 当 XX 为 单调 递 碱 序列 时 ,有 
(Tz()—zrD)) SO0j>Li. 
又 因 f 为 严格 单调 递增 函数 , 则 
(f(z0))— f(z00))) 之 0 >i 
由 (x ) 式 知 ， 
> b> (ww — ww) (f(z(7)) — fri))) 委 0， 
j=itl i=k 
故 


Df lz Pwif x2)) 
1 一 之 J 
< 


D3 D3 


灰色 缓冲 算 子 理论 及 其 应 用 


又 因 g 为 严格 单调 递增 函数 , 则 
p32 f(z(i)) 
i=k 


Duwif (x0)) 
8 <g| = 


定理 4.4 设 X 二 (zx(1),z(2),…,zx(n)) 为 非 负 的 系统 行为 数据 序列 , 且 
ZGD > 0,f 这 0,w 之 0,w'0. 权 重 向 量 分 别 为 ww 二 (ww ，…,w) 和 w 二 


(ua ,其 中 必 天 0u 天 0. 若 人 过 -区 , 即 
Wl Wl 


/ 
Ww WW i = 2 


了 为 严格 单调 递增 (或 递减 ) 函 数 ,8g 为 其 反 函 数 , 则 : 


(1) 当 X 为 单调 增长 序列 时 ,有 
g([IE f(z00)) ss )> a([TIfY Cz00)) as ); 
jk i=k 
(2) 当 X 为 单调 递减 序列 时 ,有 


sg([ 末 Ar)Jarss)<s([ 让 /rczGDD)Jaess)， 
证 明 : 先 对 函数 g 内 部 的 表达 式 两 边 取 对 数 ,再 进行 推理 ,证明 ,其 推导 过 
程 与 定理 4. 3 类 似 ,最 后 再 由 g 为 严格 单调 递增 (或 递减 函数 , 即 可 得 到 结论 . 


4.2 实例 分 析 


以 上 海 市 国际 互联 网 用 户 数 为 例 *, 验证 本 章 构造 的 弱化 缓冲 算 子 在 
GM(1,1) 模 型 预测 中 的 应 用 . 选取 该 市 2001 一 2007 年 国际 互联 网 用 户 数 ( 单 
位 :万 户 ) 作 为 原始 数据 : 

X = (310,420,432,633,803,957,1080). 
从 原始 数据 可 以 发 现 , 上海 市 上 网 用 户 数 增长 势头 很 猛 , 年 平均 增长 率 为 


4 


19. 52% ,如 此 高 的 增长 率 不 可 能 一 直 保持 下 去 . 因此 ,直接 用 原始 数据 建 模 , 预 
测 结果 令 人 难以 相信 . 经 过 分 析 , 笔 者 认为 上 海 市 二 网 用 户 数 与 Internet 的 刚 
刚 兴 起 以 及 政府 的 大 力 推广 等 因素 有 关 , 因 此 要 进行 若干 年 后 上 网 用 户 数 的 预 
测 ,必须 弱化 其 增长 趋势 ,将 上 述 政策 因素 附加 到 过 去 的 年 份 中 ,从 而 消除 前 期 
政策 因素 对 后 期 上 网 用 户 数 增加 速度 的 影响 ,使 得 模型 预测 精度 更 高 , 预测 结 
果 与 实际 情况 相符 合 . 
以 2001 一 2006 年 的 数据 作为 建 模 数据 ,以 2007 年 的 数据 作为 模拟 检验 数 

据 . 为 了 方便 , 令 f(z) = x*' ,权重 为 等 权重 ,利用 D; 和 D, 对 原始 序列 进行 组 
冲 作用 ,得 到 一 阶 缓冲 序列 : 

XD;, = (639. 21,684. 96,735. 27,809. 63,883. 70,957)， 

XD, = (548. 94,615. 39,677. 06,786. 46,876. 62,957). 
通过 计算 得 到 平均 相对 误差 和 一 步 预测 误差 ,如 表 4. 1 所 示 . 

表 4.1 弱化 前 后 模型 的 平均 相对 误差 和 一 步 预测 精度 比较 


模型 弱化 缓冲 算 子 平均 相对 误差 (%) | 一 步 预测 误差 (%) 
1 无 5.732 11.947 
2 XD, 0. 451 3. 410 
3 A 1.268 0. 162 


由 表 4.1 可 以 看 出 ,原始 序列 经 过 缓冲 算 子 D, 作用 后 ,一步 预 测 精度 最 
高 ,其 预测 模型 为 : 
(2001+1t) = 5173. 731en nz2sasx 一 4624. 7955. 
得 2007 年 上 海 市 互联 网 用 户 数 的 预测 值 为 1081. 747 万 户 ,与 实际 值 基本 
吻合 . 
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基于 单调 函数 的 新 强化 缓冲 算 子 理论 研究 


5.1 强化 缓冲 算 子 序列 


刘 思 峰 、 党 次 国 等 教授 在 其 专著 (灰色 系统 理论 及 其 应 用 ) 中 构造 了 下 列强 
化 缓冲 算 子 , 设 X = (z(1),z(2),…,z(Cz)) 为 系统 行为 数据 序列 , 令 
XD; = (z(1D)di za)di) ,i = 1,2, 
其 中 


(wr 十 … 十 wh)z Ck) 
Wwer(k) 二 十 wrx (n)” 


x (k) 二 
[zm (CA) X … X ze (n) J 
则 当 X 为 单调 增长 序列 ,单调 衰减 序列 或 振 功 序 列 时 ,，D，、D， 皆 为 强化 缓冲 算 
子 . 


Zz(k)di = 


x(k)d; 


1，2，…m， 


在 此 ,我 们 在 强化 缓冲 算 子 D, 、D 的 基础 上 ,利用 单调 函数 理论 构建 新 的 
强化 缓冲 算 子 . 
定理 5.1 设 X= (zr(1),x(2),…,Xx(n)) 为 系统 行为 数据 序列 , 目 zx(i) 之 0， 
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了 0,wi 0 .权重 向 量 志 二 《wi，…,w,), 了 为 严格 单调 递增 (或 递 碱 ) 函数 ， 
8 为 其 反 函数 . 令 

XD; = (zx(Dds,***,z(n)ds), 
其 中 


1,2,%,n, 


x(k)ds al (wrt et ww) f(x(k)) |， 


Wf (TR) Fe Ff zn)) 
则 当 X 为 单调 增长 序列 ,单调 衰减 序列 或 振 功 序列 时 ，D 为 强化 缓冲 算 子 . 
证 明 : 容 易 验 证 


/wf zn) 
zd 一 人 玉生 太一 (00， 


即 Ds 满足 缓冲 算 子 公理 一 . 
至 于 缓冲 算 子 公理 二 .公理 三 ,D; 显然 满足 . 因而 , D; 为 缓冲 算 子 . 


假设 f 为 严格 单调 递增 函数 . 
(1) 当 X 为 单调 增长 序列 时 ,有 
0<zx(k) rn), 
则 
0 一 Atz(D) < fn), 
PN or Od 
Cu + ww) fz(h)) 
Of Te Ff rn)) < fr)), 
Lo ee ee -nn 


所 以 , D; 为 强化 缓冲 算 子 . 
(2) 当 X 为 单调 衰减 序列 时 ,有 
Ze) rn) 之 0， 
则 
f(z(k)) fz) >0, 
CW tw) fzR)) wf (zh) t+ wf (rn)) > 0, 


台 灰色 缓冲 算 子 理论 及 其 应 用 加 


(we + 十 to 大 CC 
wf (TR)) + 二 wf (rn)) 宇 f(r(k)) > 0, 


Cu tt) fr)) = 
(a o> ef) i 


所 以 , D; 为 强化 缓冲 算 子 . 
(3) 当 XX 为 振荡 序列 时 , 令 
ZX(k) = mei) ,zh) = Bax), 
对 任意 iE {1,2,,…,n}), 有 
Z(R) SID, Tn) Th) Cr), rn). 
fz(k)) fr ) ,flr(n)) > 0, 
0<flrh)) fx) ,flr(n)), 


xr(kds=g 


则 
0 < wf Cr) + tf zDD) < Cor te tw ) fzh)), 
0< Cw tot ww) fh)) Swf (zhD) + tf zm)). 
Cw tet ww) f (zh)) 
Wf AR Te fr) > 7k) > 0, 
(wk 十 十 wy) f° (x(k)) 
(ed = a (mf > sz0D)) = z(h) = maxx(i); 


(wt) f(x(h)) 
Wf r+ fren) < fx)), 


0o0< 


(v0 二 十 wo)f (x(h)) x a . 
xD Gye 全 厅 j 过 sdcz0D7) 0 = mixGDi 
Maxz(i)d, 宇 r(k)d; 宇 rx(k) = max (i), 


nin Z(i)ads rh)d;s rh)= sin I(D). 


故 D; 为 强化 缓冲 算 子 . 
同 理 可 证 , 当 f 为 严格 单调 递减 函数 时 , D; 也 为 强化 缓冲 算 子 . 
定理 5.2 设 X=(z(1),x(2),…,x(n)) 为 系统 行为 数据 序列 , 且 zx(i) 之 0， 
了 >0uu 二 0. 权 重 向 量 记 一 (ww…zw) ， 太 为 严格 单调 递增 (或 递减 ) 函 数 ， 
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8 为 其 反 函 数 . 令 
XD, = (Z(1)d wzZCza)d)， 
其 中 
PCzCD) ce 
2 li (zk)) XX fo (zm) J 上 2 


则 当 X 为 单调 增长 序列 .单调 衰减 序列 或 振荡 序列 时 , D, 为 强化 缓冲 算 子 . 
证 明 : 容 易 验证 


站 f° (zr(n)) 2 
z(md, = 2 人 |= zco， 


即 D, 满足 缓冲 算 子 公理 一 . 
至 于 缓冲 算 子 公理 二 公理 三 ,D, 显然 满足 . 因而 , D, 为 缓冲 算 子 . 


假设 f 为 严格 单调 递增 函数 . 
(1) 当 X 为 单调 增长 序列 时 ,有 
0<zx(k) rn), 
则 
0< f(r) ZS frm), 
0 到 (zi) CLF Crk)) X…X fo (x0) J , 
(x(k)) 
0< £ < f(z)), 
[Pa Cz(R)) XX fo (x1)) J Se 
Cz()) 
的 sl (x(k)) ee J< ee) = a 


所 以 , D, 为 强化 缓冲 算 子 . 
《2) 当 X 为 单调 衰减 序列 时 ,有 
zk) 之 … 二 zGoD 二 0， 
则 
fz) fr) 之 0， 
f(z)) 三 [Pa (rR) X…X f(z)) a > 0, 


) 灵 色 缓冲 算 子 理 论 及 其 应 用 二 


3 二 Fz(o) >0, 
[Fa (z(D) XX fx) rT 6 


f(x(k)) 
Lpa (zk) X…X fxm)) rs 


所 以 , D, 为 强化 缓冲 算 子 . 
(3) 当 X 为 振荡 序列 时 , 令 
ZX(k) = maxz(i,z(h) 一 min r(D), 
1<isn 1Sien 
对 任意 iE {1,2,,…,n}, 有 
ZX(k) I Tn) Th) SID ,rn). 

fr(k)) f(r), fxn)) > 0, 
0<f(zh)) < fz) ,oe, fxn)), 


zx(k)d, 一 sl }> Bf rR))) = zk). 


则 
0 一 [fo Cr)) Xe Xf zn)) Tams < frCk)), 
0 二 F(z(D) 委 [Fa Cz)) X…X fr zm) ms. 


f(z) 
[Pa (zk)) XX fh (zm)) ITs 


Zz(k)di > g(f (rR))) = rR); 


f(r(k)) >0, 


< fx(h)) < fr) 
[Pa (Cz) X… X f(r)) J ef 


zh)d: < gf (zh))) = rh)s 


max zi)d, 宇 r(k)d, 二 工 (R) = Max (i), 
min xz(i)d; <r(h)d, <rh) = min zx(i). 
TS 1<in 
故 D, 为 强化 缓冲 算 子 . 
同 理 可 证 , 当 / 为 严格 单调 递减 函数 时 , D, 也 为 强化 缓冲 算 子 . 
定理 5.3 设 X= (z(1),z(2)，……z(nz)) 为 系统 行为 数据 序列 , 且 zGD 之 0， 
了 0 二 0uf> 0. 权重 向 量 分 别 为 记 一 (mm 和 am) 和 w 二 (wh ,wh)， 


第 5 章 “基于 单调 函数 的 新 强化 缓冲 算 子 理论 研究 大 


7 TEN 人 
其 中 zwm 关 0,rwn 天 0. 者 元 “ 之 z 即 


7 
WW WW i9i = 2 


J 为 严格 单调 递增 (或 递减 ) 函 数 ,8 为 其 反 函 数 , 则 ; 


(1) 当 X 为 单调 增长 序列 时 ,有 
usezdD) Dw’f: zk)) 
es <g| 气 3 
Dwif C(x)) Dwif (xr)) 
名 各 
(2) 当 X 为 单调 递减 序列 时 ,有 
Fw fi (zk) Dw fzk)) 
8| 尝 二 8| 学 2 
Dwif (zi)) Duwif ri)) 
i i=k 
证 明 : 假 设 f 为 严格 单调 递增 函数 , 则 g 也 为 严格 单调 递增 函数 , 欲 证 
wz(zGo) Swifi zk)) 
<gl 芋 ， 
Duwif (zx)) Duwif xD)) 


则 只 需 证 明 


Doif zx) Dw fi Crk)) 
Fd 冬 


Dr) ~ Di) . 
考虑 D3 Dt?) 过 Dw Df) 的 符号 . 
通过 计算 ,有 So 
PD’ Dt) 一 Dw Duif ead) 
四 四 GD 
十 … 十 zolrokf(Cz(E)) 十 rzokh f (rk 二)) + + ww,f (rn))] 


灰色 缓冲 算 子 理论 及 其 应 用 


— [ew f (zh) + wi flrlk+ D) TT wf (zm) + wf (zh)) 
十 和 … 十 zwrokrCz(R)) 十 wryf (zk 十 1)) 十 … 十 wwsf(xCn))] 


二 了 BD) ww ww) fx) — Fr)). (*%) 


jhj= 计 1 


由 


/ 
TU TUi 

> ,i = 2,3,. Nn, 
Wi Wi 


得 


i Wi » YUF 
wwXH ,wi > wxX 
Wl Wj-2 


ww WX > Xj >i 
和 


3 
则 
(ww’— wrw’) 0, > i 
(1) 当 X 为 单调 增长 序列 时 ,有 
(zx) 一 z(GD) 0,j>i 
又 因 f 为 严格 单调 递增 函数 , 则 
(f(r) — f(r) 过 0 >i 
由 (* ) 式 知 ， 


和 > ww ww) (fz) — f(r))) > 0, 


4 二 1 


ui2z(Cz(D)) > ozCzCD)) 

所 < 宇 . 

Dwf zDD) Duwif zi)) 
入 


i=k 
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又 因 g 为 严格 单调 递增 函数 , 则 
Dw fr (zh)) FF? zk) 
放生 A . 
Dwif (xi)) Dwif zi)) 


《2) 当 X 为 单调 递减 序列 时 ,有 
(xz —zD) SO0j>i 
又 因 了 为 严格 单调 递增 函数 , 则 
(fr) — fx) SO0j>i 
由 (* ) 式 知 ， 


> > ww www) fz) — f(r))) <0, 


4 一生 1 


wzce) Df (zk)) 
= > 学 . 
Dwf lr) Dwif zi)) 
er Ed 

又 因 8 为 严格 单调 递增 函数 , 则 


Dwifr (zc))) po 
可 所 兰 g| 学 
| Dwf ez)) | Pw’f zi)) 
Ed Ed 
定理 5.4 设 X= (zx(1) ,zx(2),…,z(n)) 为 系统 行为 数据 序列 , 且 zx(i) 之 0， 
>0,w; 之 0,rwf>0. 权 重 向 量 分 别 为 w= (wwww) 和 w 二 (wi ww) ， 


hn 
其 中 性, 关 0ww 关 0. 车 于- 过 -区 -， 即 
Wl Wi 
Ww Ww ii = 2 


了 为 严格 单调 递增 (或 递减 ) 函 数 , 8 为 其 反 函 数 , 则 : 
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(1) 当 X 为 单调 增长 序列 时 ,有 
_f (Crk) f(x(k) 
® [Hm SE i 
(2) 当 X 为 单调 递减 序列 时 ,有 
f(z(k)) f(z(k)) 
dH (7(0)) J B 8 a 


证 明 : 先 对 函数 g 内 部 的 表达 式 两 边 取 对 数 ,再 进行 推理 ,证 明 ,其 推导 过 
程 与 定理 5. 3 类 似 ,最 后 再 由 g 为 严格 单调 递增 (或 递减 ) 函 数 , 即 可 得 到 结论 . 


5.2 实例 分 析 


以 人 均 电力 消费 量 (单位 :千瓦 时 ) 为 例 来 验证 强化 缓冲 算 子 在 GM(1,1) 
预测 过 程 中 的 作用 . 选取 2000 一 2006 年 中 国人 均 电 力 消 费 量 作为 原始 数据 序 
列 ,2000 一 2006 年 中 国人 均 电 力 消 费 量 如 表 5. 1 所 示 . 

表 5.1 中 国人 均 电力 消费 量 (单位 :千瓦 时 ) 


| 年 份 2000 | 2001 | 2002 | 2003 | 2004 | 2005 | ze ] 
人 均 电力 消费 量 | 132.4 144.6 156. 3 173.7 190.2 216.7 249.4 
目前 ,中 国 处 在 经 济 发 展 的 工业 化 进程 中 ,经 济 的 发 展 对 电力 的 依赖 比较 
大 . 自从 2000 年 以 来 ,中 国 克服 了 亚洲 金融 危机 的 影响 ,政府 采取 了 积极 的 财 
政 政策 和 稳健 的 货币 政策 ,为 经 济 增长 注入 了 活力 ,引起 了 电力 需求 的 激增 . 因 
此 ,本 文 以 2000 一 2005 年 中 国人 均 电 力 消费 量 的 原始 数据 序列 作为 建 模 数 据 ， 
2006 年 的 数据 作为 模型 检验 数据 . 计算 得 ,人 均 电 力 消费 量 环比 增长 率 依次 为 : 


9.215% 、8. 091%、11. 132%、9. 499、13. 933% 、15. 090%, 年 平均 增长 率 为 
9.468%. 对 原始 数据 序列 进行 准 光滑 性 检验 ,计算 可 得 , 当 上 之 2003 时 ,光滑 比 


全 入- 
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分 别 为 :0. 401、0. 313、0. 272.0. 246, 均 在 (0,0. 5) 内 , 且 光 滑 比 递减 ,满足 准 光 
滑 性 条 件 . 因此 ,原始 序列 的 一 次 累加 生成 序列 具有 准 指 数 规律 中 . 但 是 明显 看 
出 ,原始 序列 前 半 部 分 增长 速度 较 慢 , 而 后 半 部 分 增长 速度 较 快 . 因此 ,最 好 对 
原始 数据 序列 进行 平滑 处 理 ,以 削弱 冲击 扰动 因素 的 干扰 ,凸显 数据 的 规律 性 . 

为 此 , 取 f(x) = g(x) 一 工 来 构造 缓冲 算 子 ,以 本 章 构造 的 缓冲 算 子 对 原 
始 数据 进行 二 阶 强化 处 理 ,建立 预测 模型 ,并 和 原始 数据 序列 直接 建 模 进行 比 
较 , 见 表 5. 2. 

表 5.2 GM(1,1) 模 型 的 预测 值 及 相对 误差 


GM(1,1) 模 型 预测 值 | 预测 相对 误差 /% 


全 (2000 十 1) = 1317. 848e* 2 一 1185. 448 236. 831 5. 040 


全 (2000 十 0) = 680. 586en 5 一 569. 765 237. 499 4.772 


(2000+ 4) = 821. 431en sa 一 706. 111 241. 210 3. 284 


由 表 5. 2 知 ,原始 数据 序列 经 过 强化 缓冲 算 子 Di 、D; 二 次 作用 后 ,所 得 到 
的 强化 缓冲 序列 都 比 原始 数据 序列 直接 建 模 的 预测 相对 误差 小 ,其 中 D; 作用 
后 得 到 的 二 阶 强化 缓冲 序列 的 预测 相对 误差 最 小 ,预测 值 为 241. 210, 比 较 逼 近 
观测 值 249. 4, 一 步 预 测 相对 误差 只 有 3. 284% , 预测 精度 较 高 . 

取 f(x) 二 式 ,g(z) 二 x 来 构造 缓冲 算 子 , 以 本 章 构造 的 缓冲 算 子 对 原始 
数据 进行 二 阶 强化 处 理 ,建立 预测 模型 ,并 和 原始 数据 序列 直接 建 模 进 行 比较 ， 
见 表 5. 3. 

表 5.3 GM(1,1) 模 型 的 预测 值 及 相对 误差 


序列 GM(1,1) 模 型 预测 值 | 预测 相对 误差 /% 


人 
X | 人 ouw+D = 9210, lleeaa — 8714. 21 236. 831 5. 040 


XD3 | 全 人 十 1) = 883. 42en ms 一 765.72 239.8 3. 84 
XD; | 全 (十 1D) = 8491. 94enas 一 399. 44 252.7 1.36 
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由 表 5. 3 知 , 原 始 数据 序列 经 过 强化 缓冲 算 子 D:、D, 二 次 作用 后 ,所 得 到 
的 强化 缓冲 序列 都 比 原始 数据 序列 直接 建 模 的 预测 相对 误差 小 ,其 中 D, 作用 
后 得 到 的 二 阶 强化 缓冲 序列 的 预测 相对 误差 最 小 ,预测 值 为 252. 7 ,比较 逼近 观 
测 值 249. 4, 一 步 预测 相对 误差 只 有 1. 36%% ,预测 精度 较 高 . 
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6.1 缓冲 算 子 的 性 质 


定理 6.1 设 X=(z(l),z(2),…,z(z)) 为 非 负 的 系统 行为 数据 序列 , 且 
ZX(i) 二 0, > 0, 为 严格 单调 递增 函数 ,8g 为 其 反 函 数 . 若 忆 是 由 工 1)，…， 
Z()，…z(D) 所 组 成 的 强化 缓冲 算 子 ,在 d 中 ,用 f(z(k)) 替换 zk)( 一 1，…z) ， 
对 所 得 到 的 表达 式 用 函数 g 去 作用 , 则 新 表达 式 e 也 为 强化 缓冲 算 子 . 

证 明 : 因 为 d 为 强化 缓冲 算 子 , 则 

Za)d 一 工 (zy 
f(zr(n)d = f(r(n)), 
则 
xr(We=g{f(r(n)d}= g(f(r(m)) = x(n), 
即 e 满 足 缓冲 算 子 公理 一 . 

至 于 缓冲 算 子 公理 二 、 公 理 三 ,e 显然 满足 . 因而 , e 为 缓冲 算 子 . 

假设 f 为 严格 单调 递增 函数 ,下 证 : 
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(1) 当 X 为 单调 增长 序列 时 ,有 
0=r(k) Zr, 


0< f(z(k)) < < flrn)). 
由 4 为 强化 缓冲 算 子 ,得 
ZC)d < zr(k), 


JCzCe))d < frk)). 
又 g 为 严格 单调 递增 函数 , 则 有 
Xx(k)e = g(f (zk))d) <gCFCzCR))) = rk), 
所 以 e 为 强化 缓冲 算 子 . 
(2) 当 X 为 单调 衰减 序列 时 ,有 
XT(k) 宇 "…* 宇 z(n) > 0, 
得 
f(x(k)) f(rn)) > 0. 
由 d 为 强化 缓冲 算 子 ,得 
rx(k)d > zr(k), 


fz(k))d > f(r(k)). 

又 g 为 严格 单调 递增 函数 , 则 有 

Xx(k)e= g(f(r(k))d) > g(f(r(k))) = rk), 

所 以 。 为 强化 缓冲 算 子 . 
(3) 当 X 为 振荡 序列 时 , 令 
ZX(k) = Max zi) ,zh) = Min zi), 
Xk) I) ,rn) rh) Sr) ,rn), 
fx(k)) f(z)) ,flrln)), 
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f lzh)) fz) ,oe, fxn)), 
f(r(k)) = maxf (x(2i)), f(r(h)) = Diaf (xz)), 
由 d 为 强化 缓冲 算 子 ,得 


maxz(i) < maxr(Dd, min x(i) > min r(Dd, 
lSien ln ISicn Iien 


故 
maxf (zx(2)) 二 maxf (zx(i))d, 
Dinf (x2)) > migf (z(i))d. 
又 g 为 严格 单调 递增 函数 , 则 有 
maxz(i)= g(rmaxf (x(i))) 
< g(maxf (zi))d) 一 max g (f(x(i))d) 
= maxr(k)e, 
im 
Mn)= Bg(minf (x(i))) 
> g(minf (x(i))q) 一 Ming (f(z(i))d) 
= minz(A)e. 
lSicn 
所 以 。 为 强化 缓冲 算 子 . 


定理 6.2 设 X= (z(1),z(2),……z(z)) 为 非 负 的 系统 行为 数据 序列 , 且 
xi) 盖 0, 太 > 0, 了 为 严格 单调 递增 函数 , g 为 其 反 函 数 . 若 d 是 由 工 (1)，…， 
ZX(k),…,r(n) 所 组 成 的 弱化 缓冲 算 子 ,在 d 中 ,用 f(x(k)) 替换 xz(k)(k 一 
1,…,n) ,对 所 得 到 的 表达 式 用 函数 g 去 作用 , 则 新 表达 式 e 也 为 弱化 缓冲 
算 子 . 

证 明 :因为 d 为 弱化 缓冲 算 子 , 则 

Xnd = x(n), 
flzn)d = f(r(m)), 
xr(me= g{f(r(n))d}= g(f(r(m)) = zr(n), 
即 e 满 足 缓 冲 算 子 公理 一 . 
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至 于 缓冲 算 子 公理 二 .公理 三 ,e 显然 满足 . 因而 , e 为 缓冲 算 子 . 
假设 f 为 严格 单调 递增 函数 ,下 证 : 
(1) 当 X 为 单调 增长 序列 时 ,有 

0<z(k) < rn), 


0<< (ze) < < frm)). 
由 d 为 弱化 缓冲 算 子 ,得 
ZI(k)d > x(k), 


fxr(k))d > f(r(k)). 
又 8 为 严格 单调 递增 函数 , 则 有 
Tk)e= g(f(r(k))d) > g(f(r(k))) = rk). 
所 以 。 为 弱化 缓冲 算 子 . 
(2) 当 XX 为 单调 衰减 序列 时 ,有 
Z(A) 之 … 二 Zr) > 0, 


fr(k) f(r) 二 0. 
由 d 为 弱化 缓冲 算 子 ,得 
zx(k)d Zz(k), 


fr(k))d < flrCk)). 
又 & 为 严格 单调 递增 函数 , 则 有 
Zz(k)e= g(f(r(k))d) < gf (zk))) = rk). 
所 以 。 为 弱化 缓冲 算 子 . 
(3) 当 义 为 振荡 序列 时 , 令 


ZX(k) = maxr(i),r(h) = min rx(i), 
leien Tsien 
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Z(A) I) ,Tn 7h) Cx), zn). 
f(z(k)) fr)) ,fr(n)), 
fz) SAFI) ,oe, fxn)), 
flr(k)) = maxf (x(D)), flrCh)) = Sinf (x7)). 
由 4 为 弱化 缓冲 算 子 ,得 


maxz (i) > maxz(Dd, min 工人 < minz (Dd, 


maxf (zx(i)) > maxf (x(i))d, 
lien li 
minf (x(D)) < minf (x(D))a. 
又 g 为 严格 单调 递增 函数 ， 
maxz(i)= B(maxf (xz(i))) 
> 8(maxf (zi))d) = Maxg (f(z(i))d) 
= maxzr (k)e, 
lSiSn 
Win zi)= &(minf (zi))) 
< g(minf (zx(i))d) = Din 8 (f(z(i))d) 
= min x(k)e. 


icn 
所 以 e 为 弱化 缓冲 算 子 . 


6.2 实例 分 析 


【实例 1】 某 企业 开发 出 一 种 新 产品 ,其 月 销售 额 如 表 6. 1 所 示 . 
表 6.1 新 产品 的 月 销售 额 (单位 :万 元 ) 


月 份 1 ol | a | |. 7 
销售 额 60.8 | 62.6 | 65.7 | 70.4 | 77.4 | 86.7 | 98.6 
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由 表 6. 1 可 以 计算 出 该 企业 新 产品 各 月 销售 额 的 增长 率 分 别 为 2. 96%， 
4.95%,7. 12% ,9.94%,12.01%,11.65%. 由 此 可 以 看 出 ,原始 数据 序列 前 半 
部 分 的 增长 速度 明显 慢 于 后 半 部 分 的 增长 速度 . 通过 统计 分 析 发 现 ,该 企业 从 3 
月 份 开 始 采取 了 一 系列 促销 手段 ,所 以 从 4 月 份 开始 该 新 产品 的 月 销售 额 增长 


率 比 1 


一 3 月 份 有 了 明显 的 增长 . 因此 ,为 了 提高 8 月 份 及 其 以 后 月 份 的 销售 额 


预测 精度 ,需要 利用 强化 缓冲 算 子 对 原始 数据 序列 进行 强化 作用 ,以 消除 促销 


手段 对 


原始 数据 序列 的 冲击 扰动 . 


下 面 用 前 6 个 月 的 销售 额 为 建 模 数 据 ,并 预测 7 月 份 的 销售 额 . 为 了 提高 


GM(1 


,1) 模 型 的 预测 精度 ,必须 对 原始 数据 序列 进行 强化 处 理 , 用 强化 缓冲 算 


子 对 原始 数据 序列 进行 强化 处 理 ,使 得 预测 结果 与 实际 情况 更 为 接近 . 
利用 强化 算 子 ( 这 里 取 f(z) 二 懂 , 则 其 反 函 数 g(x) = Vx ) 对 原始 数据 序 
列 进 行 一 阶 强化 处 理 , 得 到 强化 数据 序列 为 


XI = (51.95,53. 64,57. 20,63. 18,72. 90,87. 6). 
表 6.2 两 种 情况 GM(1,1) 模 型 


序列 预测 值 预测 相对 误差 
直接 建 模 92. 62 4.3% 
XE 96. 66 一 0.14% 
-| | 


表 6.3 基于 强化 缓冲 算 子 的 建 模 误差 检验 表 


序 号 实际 数据 模拟 数据 相对 误差 
1 53. 64 50. 56 ”5.58% 
2 57. 20 57. 64 0.76% 
| 3 63. 18 65. 60 3.82% 
4 72.90 74.64 2. 39% 
5 87. 60 84. 94 —3.04% 
6 96. 80 96. 66 一 0. 14% | 
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月 份 的 预测 销售 额 为 96. 66, 与 实际 的 98. 6 相差 不 多 , 相对 误差 仅 为 


一 0. 14%, 比 直接 建 模 的 预测 精度 高 . 


= 
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【实例 2〗 1997 年 以 后 , 随 着 亚洲 金融 危机 对 我 国 的 影响 逐步 减弱 ,以 及 
国家 扩张 性 财政 政策 的 有 效 实施 和 消费 结构 的 升级 ,2002 年 我 国 GDP 已 突破 
10 万 亿 元 大 关 , 经 济 发 展 进入 新 一 轮 的 增长 周期 . 从 2003 年 开始 ,能 源 消耗 量 
急剧 增长 ,1998 一 2006 年 我 国 能 源 消耗 总 量 的 原始 数据 如 表 6. 4 所 示 . 

表 6.4 1998 一 2006 年 我 国 能 源 消耗 总 量 (单位 : 亿 吨 ) 


13. 22 | 13. 38 | 13. 86 | 14. 32 | 15. 18 | 17. 50 | 20. 32 | 22. 47 | 24. 63 


由 表 6. 4 可 计算 出 ,1999 一 2006 年 我 国 能 源 消耗 增长 速度 分 别 为 :1. 2%， 
3.5%,3.4%,6.1%,15.3%,16.1%,10.6%,9.6%. 2003 年 以 前 的 增长 速度 显 
著 慢 于 后 半 部 分 的 增长 速度 . 如 果 直 接 基于 原始 数据 进行 建 模 , 则 2003 一 2006 
年 的 平均 预测 误差 竟 高 达 一 80% 以 上 ,这 样 的 结果 是 不 能 接受 的 . 为 了 及 时 准 
确 地 把 握 能 源 消耗 的 发 展 趋势 ,对 能 源 需求 进行 科学 合理 的 预测 , 必须 对 缓慢 
增长 的 数据 加 以 处 理 , 使 其 符合 2002 年 后 的 发 展 趋势 ,在 此 基础 上 进行 合理 的 
预测 . 

下 面 以 前 8 年 的 数据 作为 建 模 数据 ,并 预测 2006 年 的 消耗 量 . 为 了 提高 
GM(1,1) 模 型 的 预测 精度 ,必须 对 原始 数据 序列 进行 强化 处 理 , 用 强化 缓冲 算 
子 对 原始 数据 序列 进行 强化 处 理 , 使 得 预测 结果 与 实际 情况 更 为 接近 . 

利用 强化 算 子 ( 这 里 取 f(x) 二 式 , 则 其 反 函数 g(x) = VX ) 对 原始 数据 序 
列 进行 一 阶 强化 处 理 , 得 到 强化 数据 序列 为 

Xi = (10. 53,10. 51,10. 94,11.26,12.08,15. 16 ,19. 27,22.47). 
表 6.5 两 种 情况 GM(1,1) 模 型 


序列 预测 什 | 预测 相对 误差 
直接 建 模 24.09 | 一 2.2% 
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表 6.6 基于 强化 缓冲 算 子 的 建 模 误差 检验 表 


模拟 数据 相对 误差 。 ] 
8.76 —16. 68% 

时 10. 94 10. 17 —7.02% | 
3 11.26 11.82 4.94% 
13.63% 


5.19% 
—3. 88% 
—4.25% 


2006 年 的 预测 消耗 量 为 24. 99, 与 实际 的 24. 63 相差 不 多 ,相对 误差 仅 为 
1. 5%, 比 直接 建 模 的 预测 精度 高. 
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构造 及 其 应 用 


7.1 一 类 新 的 强化 缓冲 算 子 的 构造 


定理 7.1 设 X= (x(1),zx(2),…,z(n)) 为 系统 行为 数据 序列 ,z(A) 之 0， 


和 一 1,2,…,n ,其 缓冲 序列 为 XD, 一 (z(1)di ,rz(2)d yz(a)di) ,其 中 


ze = 1,2,0 ns 


x(n) 


(di = x(k) [ 


则 当 XX 为 单调 增长 序列 .单调 衰减 序列 或 振荡 序列 时 , D, 都 为 强化 缓冲 算 子 . 


证 明 : 容 易 验 证 , D, 满足 缓冲 算 子 的 三 个 公理 ,所 以 Di 为 缓冲 算 子 . 下 证 : 


Di 为 强化 缓冲 算 子 . 
(1) 当 X 为 单调 增长 序列 时 , 则 


Es eg 
zd —z(h) = zk) [EY] — zk) 


一 xD[[ 包 | —1]<o, 
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因此 
zd SrA), 
由 定理 1.1~1.3 知 ,缓冲 算 子 Di 为 强化 缓冲 算 子 . 
(2) 当 X 为 单调 衰减 序列 时 ,有 
(A) 


Ea 
zd — z= zh [EE] ”一 区 


.= zB TS 一 1] 宕 0， 


x(n) 
因此 
x(k)di 之 工 (R)， 
由 定理 1. 1 一 1. 3 知 ,缓冲 算 子 D, 为 强化 缓冲 算 子 . 
(3) 当 X 为 振荡 序列 时 , 令 
x(a) = max{z(k)|k= 1,2,.%,n), 
Xx(B) = min{x(k)|k= 1,2,%,n)s 


zx(a) 2 
Xx(n) 


max x (Dd 宇 z(a)di 一 Za) [ 宇 x(a) = maxz(i), 


minz (Dd < z(Pd 一 (有 [8]™ < (8) = min zi). 
由 定理 1. 1 一 1, 3 知 ,缓冲 算 子 Di, 为 强化 缓冲 算 子 . 
定理 7.2 设 X= (z(1),z(2),…,z(z)) 为 系统 行为 数据 序列 ,z(A) 之 0， 
= 1,2,…,n ,其 缓冲 序列 为 XD。 = (z(l1)dosz(2)ds，yz(z)dn)》，mma 为 自 
然 数 ,其 中 


zh)d, 二 工 (有 E23 ss 
则 当 XX 为 单调 增长 序列 ,单调 衰减 序列 或 振荡 序列 时 ，D,, 都 为 强化 缓冲 


算 子 . 
证 明 : 容 易 验证 , D, 满足 缓冲 算 子 的 三 个 公理 ,所 以 D。 为 缓冲 算 子 .下 证 ， 
了 , 为 强化 缓冲 算 子 . 
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(1) 当 六 为 单调 增长 序列 时 ,有 


be 
zk)dn — x(k)= x(k) (Ee 


因此 
IT(k)d» < x(k), 
由 定理 1. 1 一 1. 3 知 ,缓冲 算 子 D, 为 强化 缓冲 算 子 . 
(2) 当 X 为 单调 衰减 序列 时 ,有 


A 
zd —z(h)= zk) (EB) 一 


| Zz TE 
=z([ (E07)) 1]>0, 


因此 
I(k)dn > rk), 
由 定理 1. 1 一 1. 3 知 ,缓冲 算 子 D, 为 强化 缓冲 算 子 . 
(3) 当 X 为 振荡 序列 时 , 令 
xz(a) = max{z(k)|k = 1,2,.…,n}, 
Zz(B) = min{z(k) Ik = 1,2,..,n}s 


x(W dn = x(a) (Ee 


EC 过 za)， 


a 
z(pd = 7(B) (EB) <zp, 
因此 
max (Ddn 宇 z(a)d», 三 Za) = Max), 


min z(i)du < zr(Bd, <zx(P = min ri). 
1<i<n lSicn 


由 定理 1. 1 一 1. 3 知 ,缓冲 算 子 D, 为 强化 缓冲 算 子 . 


天 次 色 缓冲 算 子 理 i 及 其 应 用 


推论 1 当 m 二 nn 一 十 1 时 ,得 到 文献 [22] 定 义 的 强化 缓冲 算 子 , 即 
XD = (z(Dd,z(2 d,s zmd), 
其 中 
zd = x(k) LK, = 1,2,..,n, 
x(n) 


则 当 六 为 单调 增长 序列 .单调 衰减 序列 或 振荡 序列 时 , D 皆 为 强化 缓冲 算 子 . 
定理 7.3 设 X= (z(1),z(2),…,z(z)) 为 系统 行为 数据 序列 ,z(A) 之 0， 
k 二 1,2,…,n, 其 缓冲 序列 为 XD; 二 (z(1)ds,z(2)ds ，…zGa)adz) ,其 中 


Zz(k) 工 (站 吕 
(ds 5 (Gy ls,2 eon, 


则 当 X 为 单调 增长 序列 或 单调 衰减 序列 时 , D; 皆 为 强化 缓冲 算 子 . 

证 明 : 容 易 验 证 , D, 满足 缓冲 算 子 的 三 个 公理 ,所 以 D; 为 缓冲 算 子 . 下 证 ， 
D; 为 强化 缓冲 算 子 . 

(1) 当 X 为 单调 增长 序列 时 ,有 


xz Ww /rx VRHT 
x(k)ds — x(k) 二 (Be x(k) 


一 _z(k) (8) 可 1] 


将 Xx(n) 
+[ (2D) -1 *…+0 1D)<o, 
因此 
x(k)d; < zx(k), 
故 D; 为 强化 缓冲 算 子 . 
(2) 当 X 为 单调 衰减 序列 时 ， 


一 z(D) Vi 

一 上 十 2 (8) A) 
二 
= 


le, 读 虽 


+[ (2 路 了 J)>0, 


x(k)ds — zx(k) 
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TCR)dz x(k), 


因此 


故 D: 为 强化 缓冲 算 子 . 
定理 7.4 设 X= (x(1),z(2),…,z(n)) 为 系统 行为 数据 序列 ,z(e) 之 0， 
上 二 1,2,…,n ,其 缓冲 序列 为 XD; = (z(1)ads,z(2)ds ,za)ds), 其 中 


2z(k) x(D) Tn 
(Dd = TD (EK) = ln, 


则 当 X 为 单调 增长 序列 或 单调 衰减 序列 时 , D; 皆 为 强化 缓冲 算 子 . 

证 阴 : 容 易 验证 , D; 满足 缓冲 算 子 的 三 个 公理 ,所 以 D; 为 缓冲 算 子 .下 证 : 
DD; 为 强化 缓冲 算 子 . 

《1) 当 X 为 单调 增长 序列 时 ,有 


员 2x(k) ;1ZD mn 
zd 一 xD 一 ir 2 (KH) zh 


2x(k) 人 [(2@ Hn 
(nn 十 k)(n—k 十 1) x(n) 


) 
+ (k++ D[ (ED) 1+…+z0-D<o 
| ] <o. 


) 


因此 
zk)ds < zh), 
故 D 为 强化 缓冲 算 子 . 
〈2) 当 X 为 单调 衰减 序列 时 ,有 
2zx(k) /zDD) Yn 
zd rd TT (S63) 全 名 
2x(k) Zz(k) 
min RTD\ Cte) 1] 
zk+D\ 志 1 
+ (+d[( ZX(n) ) 1 +all 1)>0. 
因此 


Zz(k)d; > zx(k), 
故 D; 为 强化 缓冲 算 子 . 
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定理 7.5 设 XX 二 (x(]),z(2),…,r(n)) 为 系统 行为 数据 序列 , x(k) 之 0 ， 
k= 二 1,2,…,n, 其 缓冲 序列 为 XD, ==(z(1)4d,,z(2)d,,…,zx(n)d,) ,其 中 


xd = zk) [ 立 (Ee ji Vs 
1 


Zz(n) 
则 当 X 为 单调 增长 序列 或 单调 衰减 序列 时 ,Di 皆 为 强化 缓冲 算 子 . 
证 明 : 容 易 验 证 , D, 满足 缓冲 算 子 的 三 个 公理 ,所 以 D, 为 缓冲 算 子 . 下 证 ， 
D, 为 强化 缓冲 算 子 . 
(1) 当 X 为 单调 增长 序列 时 ,有 


za=zi [TH (E22) 
jk 


ZI(n) 


0 [EB)™ (ED zw, 


因此 
ZC)d < zk), 
故 D; 为 强化 缓冲 算 子 . 
〈2) 当 X 为 单调 衰减 序列 时 ,有 


x(k)di= zx(k) [ 立 ED Te 


x(n) 


OB) (SED) a0, 


因此 
Zz(k) ds > x(k), 
故 D, 为 强化 缓冲 算 子 . 

从 以 上 的 讨论 可 知 , 由 于 强化 组 神算 子 必 须 满足 不 动 点 公理 , 即 z(zp)d 一 z0D). 
因此 , 当 强 化 缓冲 算 子 作用 于 单调 增长 序列 时 ,数据 萎缩 , 即 强化 缓冲 序列 的 增 
长 速度 比 原始 序列 的 增长 速度 减缓 ;而 当 其 作用 于 单调 衰减 序列 时 ,数据 膨胀 ， 
即 强化 缓 串 序列 的 衰减 速度 比 原始 序列 衰减 速度 减缓 . 因此 , 当 原始 序列 的 前 
半 部 分 增长 (衰减 ) 速 度 较 缓 ,后 半 部 分 增长 (衰减 ) 速 度 较 快 时 ,可 以 利用 本 文 
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所 构造 的 强化 缓冲 算 子 对 其 作用 ,将 使 序列 变 得 比较 平缓 ,并 且 考 虑 了 “新 信息 
优先 ”的 原则 ,能 够 有 效 地 消除 冲击 扰动 对 系统 数据 序列 造成 的 “失真 ”现象 , 因 
而 能 够 提高 模型 的 预测 精度 . 


7.2 实例 分 析 


以 人 均 电 力 消费 量 (单位 :千瓦 时 ) 为 例 来 验证 本 章 所 述 强化 缓冲 算 子 在 
GM(1,1) 预 测 过 程 中 的 作用 . 选取 2000 一 2006 年 中 国人 均 电力 消费 量 作为 原 
始 数据 序列 ,如 表 7. 1 所 示 . 

表 7.1 中 国人 均 电力 消费 量 (单位 :千瓦 时 ) 


年 份 | ao | za | 2002 2003 2004 | 2005 | 2006 

人 约 电 力 消费 量 | 132.4 | 144.6 | 156.3 | 173.7 | 190.2 | 216.7 | 249.4 
目前 ,中 国 处 在 经 济 发 展 的 工业 化 进程 中 ,经 济 的 发 展 对 电力 的 依赖 比较 
大 , 自从 2000 年 以 来 ,中 国 克服 了 亚洲 金融 危机 的 影响 ,政府 采取 了 积极 的 财 
政 政 策 和 稳健 的 货币 政策 ,为 经 济 增长 注 人 了 活力 ,引起 了 电力 需求 的 激增 . 因 
此 ,本 节 以 2000 一 2005 年 中 国人 均 电 力 消费 量 的 原始 数据 序列 作为 建 模 数 据 ， 
2006 年 的 数据 作为 模型 检验 数据 . 计算 可 得 ,人 均 电力 消费 量 环比 增长 率 依次 
为 :9.215%、8.091%、11.132%、9. 499、13. 933% .15. 090% ,年 平均 增长 率 为 
9. 468%%. 对 原始 数据 序列 进行 准 光滑 性 检验 ,计算 可 得 , 当 上 之 2003 时 ,光滑 比 
分 别 为 :0. 401.0. 313、0. 272、0. 246 , 均 在 (0,0. 5) 内 , 且 光 滑 比 递减 ,满足 准 光 
滑 性 条 件 . 因此 ,原始 序列 的 一 次 累加 生成 序列 具有 准 指数 规律 1. 但 是 明显 看 
出 ,原始 序列 前 半 部 分 增长 速度 有 点 慢 ,后 半 部 分 增长 速度 较 快 . 因此 ,最 好 对 


原始 数据 序列 进行 平滑 处 理 ,以 削弱 冲击 扰动 因素 的 干扰 ,凸显 数据 的 规律 性 . 
为 此 ,以 本 章 构 造 的 缓冲 算 子 对 原始 数据 进行 二 阶 强化 处 理 ,建立 预测 


-DD 
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模型 ,如 表 7. 2 所 示 . 
表 7.2 四 种 情况 GM(1,1) 模 型 


GM(1,1) 模 型 预测 值 | 预测 相对 误差/% 
全 (2000 十 1) = 1317. 848en oz 一 1185. 448 236. 831 | 5.040 
全 (2000 十 上 = 680. 586e* xz — 569. 765 237. 499 | 4.772 
全 (2000 十 六 = 821. 431en mm 一 706. 111 241. 210 | 3.284 
人 (2000 十 z) = 876. 605er sz 一 758. 408 240. 118 | 3.722 
全 (2000 十 昌 = 817. 537e te 一 702. 348 241. 185 3.294 


由 表 7. 2 可 知 ,原始 数据 序列 经 过 强化 缓冲 算 子 D; 、D, 、D, 和 D, 作用 

~ 后 ,所 得 到 的 强化 缓冲 序列 都 比 原始 数据 序列 直接 建 模 的 预测 相对 误差 小 ,其 

中 经 Ds 作用 后 得 到 的 二 阶 强化 缓冲 序列 的 预测 相对 误差 最 小 , 预测 值 为 

241. 210, 比较 通 近 观测 值 249. 4, 一 步 预测 相对 误差 只 有 3. 284%, 即 预测 精度 

最 高 . 原始 数据 序列 X 经 过 三 阶 强化 缓冲 算 子 D; 作用 后 预测 值 为 240. 118, 预 

测 误差 为 3. 722% ;原始 数据 序列 X 经 过 四 阶 强化 缓冲 算 子 D, 作用 后 预测 值 
241. 185 ,预测 误差 为 3. 294%. 
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基于 新 信息 优先 的 弱化 缓冲 算 子 的 
构造 及 其 应 用 


8.1 一 类 新 的 弱化 缓冲 算 子 的 构造 


定理 8.1 设 X= (zx(1),x(2),…,x(n)) 为 系统 行为 数据 序列 ,z(&) 之 0， 
& 二 1,2,…,n ,其 缓冲 序列 为 XD, 一 (z(1)d ,z(2)d ,za)di) ,其 中 


= 
ER) ,k= 1,2，…m， 


CR)cd = x(k) [a 


则 当 为 单调 增长 序列 单调 衰 减 序列 或 振荡 序列 时 ,D, 皆 为 弱化 缓冲 算 子 . 
证 明 : 容 易 验证 , D, 满足 缓冲 算 子 的 三 个 公理 ,所 以 Di 为 缓冲 算 子 . 下 证 : 


DD 为 弱化 缓冲 算 子 . 
(1) 当 XX 为 单调 增长 序列 时 ,有 
zh di — zk)= zk) (Ee jy —z(k) 


-wo[ EB) -ie 
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因此 
zx(k)di > x(k), 
由 定理 1. 1 一 1. 3 知 ,缓冲 算 子 Di, 为 弱化 缓冲 算 子 . 
(2) 当 X 为 单调 衰减 序列 时 ,有 


zd —z(k)= zh)( 


zx(n) )™ 


Zk =x(h) 


因此 
I(R di < zk), 
由 定理 1. 1 一 1. 3 知 ,缓冲 算 子 D, 为 弱化 缓冲 算 子 . 
(3) 当 X 为 振荡 序列 时 , 设 
x(a) = max{I(k)|k= 1,2,..,n}), 
x(B) = min{z(k)|k = 1,2,.,n}. 
对 任意 的 i€ {1,2,…,n), 有 
zxDadi = 7 (DE (x) < x(a), 
xz(Dadi = x (DE (zm > zB). 
因此 
max {zk)d }< Max {zk)} ， 
mi {zh d}> min {7(k)}, 
由 定理 1. 1 一 1. 3 知 ,缓冲 算 子 D, 为 弱化 缓冲 算 子 . 
定理 8.2 设 X= (zx(1),z(2),…,zx(n)) 为 系统 行为 数据 序列 ,z(A) 盖 0， 
上 二 1,2,…,n ,其 缓冲 序列 为 XD, 二 《zx(1)4d;,z(2)ds,…,z(n)qds), 其 中 
EE 


4 一 二 打 凡 :0 (开放 


则 当 X 为 单调 增长 序列 ,单调 训 减 序列 或 振荡 序列 时 ，D: 皆 为 弱化 缓冲 算 子 . 
证 明 : 容 易 验证 , D; 满足 缓冲 算 子 的 三 个 公理 ,所 以 D; 为 缓冲 算 子 . 下 证 : 
DD; 为 弱化 缓冲 算 子 . 


人 == 


一 1, 2 
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(1) 当 X 为 单调 增长 序列 时 ,有 


ZO NT 
x(k)ds — zx(k) EE D (a8) zx(k) 


i=k 


[ze ye 一 z(b | 二 … 十 [zGnD —z(&)] 
n 一 k 十 1 


之 0， 


因此 
Z(A)da > x(k), 
故 D: 为 弱化 缓冲 算 子 . 
(2) 当 六 为 单调 衰减 序列 时 ,有 


zn) \ 
Zz(k)ds — x(k)= = EE (a —z(k) 


[zw (2 zk) |+ + [rn) — zk)] 
和 nk+1 0s 


因此 
x(k)d; < rlk), 
故 D; 为 弱化 缓冲 算 子 . 
《3) 当 X 为 振 功 序 列 时 , 令 
Z(a) = max{z(k)|k = 1,2,.,.n}, 
Zz(B) = min{z(k)|k = 1,2,..,n}. 
对 任意 的 yE {1,2,…,n), 有 
Hh 


zi — x(a) 2 z(0) (EP) — zo) 


[< (E28) 一 zx(o) | + [zn) —z(e)] 
x 二 


[Cz zn) 一 zao)] 十 … 十 [CzCooD 一 zx(@)] < 0， 
其 -一季 路 玉 
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即 
Z(7y)ads ra), 
因此 
Max{r(k)d}< Max {zk)}. 
同 理 可 证 
xz(y)ds > x(p), 
则 ， 


min (x(k) de} Din, (zk)}. 
故 当 X 为 振荡 序列 时 , D; 为 弱化 缓冲 算 子 . 
定理 8.3 设 X= (z(1),z(2),…,z(z)) 为 系统 行为 数据 序列 ,zx(k) 汪 0， 
二 1,2,…,n，Di 和 DD; 为 弱化 缓冲 算 子 ,如 定理 8. 1 和 定理 8. 2 所 述 , 则 
(1) 当 X 为 单调 增长 序列 时 ， 
(di Sr(k ds VERE {1,2 sn)s 
(2) 当 X 为 单调 衰减 序列 时 ， 
(di x(k dis, VkE {1,2,0,n}. 
证 明 :(1) 当 X 为 单调 增长 序列 时 ,对 Yk E {1,2,…,n 一 1} ,有 


zn)_\™ 
Zz(k+ Dd zk+D (sD) 
xT(k)d! x(n) VT 
Ck (230%) 


工 (n) | 融 
Rh (a 5) 


zw (EE)” 


a (M2) A 


即 {zxCh)di|k 二 1,2,…,n) 为 递增 序列 . 


te : 1 


人 


& 


zh) (EO) 二 zk+D (HE) + zn) 
一 此 十 1 
zx Vi Zn VRHTI 
zk) (FE) 十 … 十 CD (FE) 
wi n—k 二 1 
en 
= zk (2D) 一 Cd， 


即 
ZJ)di < zk ds. 
(2) 同 理 可 证 , 当 X 为 单调 衰减 序列 时 ,有 
x(k di > x(k) de. 
定理 8.4 设 X= (xz(1),z(2),…,z(n)) 为 系统 行为 数据 序列 ,z(A) 之 0， 
三 1,2,…,n ,缓冲 序列 为 XD; 二 (TC1Dd;,7(2)ds,…,z(n)qd;), 其 中 
(hd = TEETD 2 (6 ee 
则 当 X 为 单调 增长 序列 .单调 衰减 序列 或 振荡 序列 时 ,D， 皆 为 弱化 缓冲 算 子 
证 阴 : 容 易 验证 , D; 满足 缓冲 算 子 的 三 个 公理 ,所 以 Ds 为 缓冲 算 子 . 下 证 ， 
D; 为 弱化 缓冲 算 子 . 


,一 127 


(1) 当 X 为 单调 增长 序列 时 ,有 
Ph 
2》Viz(D (Te 
(Ada — zx(h)= 一 二 (a®) Zz(k) 


(n+R)n—k+1) 
rT nth)(n—kt Dr(k) 
2( Di ( ) aa 了 } 


zi) 
(nt k(n—k 二 1) 
2{ 丰 Zz(W Vi : 
z(k) (2 ) (6] 十 … 十 xz[z(m) z(]} 


(十 妇 ( 一 上 十 1 0 
因此 ， 

x(k) ds > zk), 
故 D; 为 弱化 缓冲 算 子 . 
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(2) 同 理 可 证 , 当 X 为 单调 衰减 序列 时 ,有 
zx(k)ds < r(k), 
故 Ds 为 弱化 缓冲 算 子 . 
《3) 当 X 为 振荡 序列 时 , 令 
za) = max{z(k) | 人 一 1,2，…m2)， 
Z(B) = min{7Y(k) | k= 1,2,..…,n}, 
对 任意 的 YE {1,2,…,n)}, 有 


/ZO Tn 
2 (8) 
(2 十 (一 7 十 1) 


ZX(Wd; — x(a) Zz(a) 


2(7z( (2) 一 za) + +zln) 一 z(g] 
加 xy, re) Q. (9 
(nn—7+) 


20[ (zh (xln)m z(t + lr) zx 0 
+n—7+1D) A 


得 
(Vd: < rla). 
同 理 可 证 
z(W ds > z(B). 
故 当 义 为 振荡 序列 时 , D; 也 为 弱化 缓冲 算 子 . 
定理 8.5 设 X= (zx(1),z(2),…,x(n)) 为 系统 行为 数据 序列 ,x(k) 二 0， 
三 1,2,…,n，D, 和 了; 为 弱化 缓冲 算 子 ,如 定理 8. 1 和 定理 8. 4 所 述 , 则 
(1) 当 为 单调 增长 序列 时 ， 
ZE)di < zr(k)ds, VEE {1,2 72)}; 
(2) 当 六 为 单调 衰减 序列 时 ， 
I(k)di 之 TOR)da YE {1l,2,.%,n). 
证 明 :(1) 当 X 为 单调 增长 序列 时 ,对 Vk E {1,2,…,n 一 1} ,有 
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工 (n)  \ 去 
+Dd zt+D (了 D) 
x(k)d za A 
Tk) (38%) 
本 
zk+D (2D) (2 DD) a 
I zk) 2 
(hk) (5 
Zz(k) 
所 以 ,，{z(k)di|k 二 1,2,…,n) 为 递增 序列 . 
2 a /zn) \H 
zd = tT mT (36) 
2 
2[ kz (EP) + Dt (= Ry) 十 … 十 mx] 
TAO 一 ATD 
Zz(n) \ TH xz) Zn) TH 
2 人 (2%) + Ck+ Dz ( Eg ) 十 十 nz(k) (5) ] 
和 (n 十 k)(n 一 k 十 1) 
zr Vi 
一 (到 的) 
= x(k)di, 
即 


x(k)di < r(k)d;. 

(2) 同 理 可 证 , 当 X 为 单调 衰减 序列 时 ,有 

x(k)di > zr(k)d:. 

从 以 上 讨论 可 知 ,由 于 弱化 缓冲 算 子 必须 满足 不 动 点 公理 , 即 x(n)4d 一 x(n) . 
因此 , 当 弱 化 缓冲 算 子 作用 于 单调 增长 序列 时 ,数据 膨胀 , 即 弱化 缓冲 序列 的 增 
长 速度 比 原始 序列 的 增长 速度 减缓 ;而 当 其 作用 于 单调 衰减 序列 时 ,数据 萎缩 ， 
即 弱化 缓冲 序列 的 衰减 速度 比 原始 序列 衰减 速度 减缓 . 因此 , 当 原 始 序列 的 前 
半 部 分 增长 (衰减 ) 速 度 较 快 ,后 半 部 分 增长 (衰减 ) 速 度 较 慢 时 ,可 以 利用 本 文 
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所 构造 的 弱化 缓冲 算 子 对 其 作用 ,使 序列 变 得 比较 平缓 ,并 且 考虑 了 "新 信息 优 
先 "的 原则 ,能 够 有 效 地 消除 冲击 扰动 对 系统 数据 序列 造成 的 “失真 "现象 ,因而 
能 够 提高 模型 的 模拟 精度 和 预测 精度 . 


8.2 实例 分 析 


以 中 国 城镇 登记 失业 人 数 (单位 :万 人 ) 为 例 来 验证 本 章 所 构造 的 弱化 缓冲 
算 子 在 GM(1,1) 预 测 过 程 中 的 作用 . 选取 2000~2006 年 中 国 城镇 登记 失业 人 
数 作为 原始 数据 序列 ,如 表 8. 1 所 示 . 
表 8. 1 ”中国 城 镇 登记 失业 人 数 (单位 :万 人 ) 


年 份 2000 2001 2002 2003 2004 2005 2006 


失业 人 数 595 681 | .770 800 827 839 847 


以 2000 一 2005 年 的 数据 作为 建 模 数据 ,2006 年 的 数据 作为 模型 检验 数据 . 
计算 可 得 ,城镇 登记 失业 增长 率 分 别 为 14. 454% .13. 069% .3. 896% .3. 375%、 
1.451% .0. 954%% ,显然 前 半 部 分 增长 速度 较 快 ,后 半 部 分 增长 速度 较 慢 , 如 果 
用 此 数据 直接 建 模 、 预 测 , 不 可 取 . 经 分 析 , 笔 者 发 现 原因 有 二 :一 方面 ,20 世纪 
90 年 代 中 后 期 到 21 世纪 初 ,由 于 国有 企业 改革 ,造成 了 许多 工人 下 岗 ; 另 一 方 
面 ,由 于 大 中 专 院 校 扩大 招生 规模 ,为 社会 培养 了 许多 的 大 学 生 , 因 此 ,就 业 压 
力 增 大 ,失业 人 数 也 大 大 增加 . 后 来 ,由 于 中 央 政 府 和 地 方 政 府 陆 续 出 台 了 许多 
促进 下 岗 职 工 再 就 业 和 扶持 大 学 生 就 业 的 政策 ,减缓 了 就 业 压 力 ,表现 为 城镇 
登记 失业 人 数 增长 减缓. 为 了 消除 原始 数据 序列 受到 冲击 扰动 因素 的 影响 ,用 
缓冲 算 子 对 其 作用 . 

以 本 章 所 构造 的 缓冲 算 子 对 原始 数据 进行 弱化 处 理 ,分 别 建立 预测 模型 
如 下 : 


OO 
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无 缓冲 算 子 作用 ,由 原始 数据 序列 直接 建立 GM(1,1) 模 型 : 
全 (2000 十 #) 一 14912. 796e*%7 一 14317. 796. 

经 缓冲 算 子 D, 作用 后 得 到 弱化 缓冲 序列 XD，, 建 立 GM(1,1) 模 型 : 
全 (2000 十 划 = 19308. 68len oa 一 18678. 607. 

经 缓冲 算 子 D: 作用 后 得 到 弱化 缓冲 序列 XD ,建立 GM(1,1) 模 型 : 
全 (2000 十 昌 = 63785. 612en ol 一 63017. 018. 

经 缓冲 算 子 D, 作用 后 得 到 弱化 缓冲 序列 XD;， ,建立 GM(1,1) 模 型 : 
从 (2000 十 上 = 98293. 850en oog 一 97490. 980. 

经 无 缓冲 算 子 和 缓冲 算 子 D 、D, 、D; 分 别 作用 后 的 缓冲 序列 所 建立 的 

GM(1,1) 模 型 的 平均 相对 误差 和 预测 值 见 表 8. 2. 
表 8. 2 ”四 种 情况 所 得 的 平均 相对 误差 和 预测 值 


平均 相对 误差 (% ) 
2.64 


2.35 


0.58 
0.26 


由 表 8. 2 知 ,原始 数据 序列 经 过 弱化 缓冲 算 子 D 、D: 和 D; 作用 后 ,平均 
相对 误差 都 比 原始 序列 直接 建 模 的 平均 相对 误差 小 ,其 中 D; 作用 后 得 到 的 弱 
化 缓冲 序列 的 平均 相对 误差 最 小 ,预测 值 为 848, 比 较 逼 近 观 测 值 847. 
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基于 有 理 插值 公式 的 GM (1,1) 模 型 
背景 值 的 构造 及 其 应 用 


9.1 GM(1,1) 动 态 预测 模型 的 建 模 机 理 


灰色 系统 理论 具有 所 需 样本 数据 少 ,不 需要 计算 统计 特征 量 等 优点 . 因此 ， 
自 1982 年 提出 以 来 得 到 了 研究 人 员 的 重视 ,已 经 在 许多 领域 ,特别 是 在 显著 不 
确定 性 和 缺乏 数据 信息 的 领域 得 到 了 成 功 应 用 . 利用 GM(1,1) 模 型 进行 预测 
虽然 有 许多 成 功 的 实例 ,但 同时 也 存在 一 些 预测 误差 过 大 的 情况 ,反映 了 
GM(1,1) 模 型 的 实用 性 有 待 提 高 . 因此 ,对 GM(1,1) 模 型 进行 深入 地 研究 , 找 
出 影响 GM(1,1) 模 型 精度 和 适应 性 的 因素 ,具有 非常 重要 的 理论 价值 和 实际 
意义 . 文献 [3] 用 实验 的 方法 分 析 了 GM(1,1) 模 型 的 误差 特性 ,文献 [4] 提 出 了 
GM(1,1) 模 型 中 的 背景 值 构造 方法 是 影响 其 精度 和 适应 性 的 关键 因素 . 基于 十 
老 的 连 分 式 理论 ,笔者 给 出 一 种 新 的 有 理 插值 ,同时 通过 仿真 例子 说 明 本 文 所 
提出 的 方法 的 有 效 性 . 

GM(1,1) 模 型 是 最 常用 的 一 种 灰色 动态 预测 模型 ,该 模型 由 一 个 单 变量 的 
一 阶 微分 方程 构成 . 其 建 模 过 程 如 下 : 


全 
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设 原 始 非 负 数据 序列 为 
Xeo = {zo (1) xz (za)}， (1) 
其 中 xm >>0,i= 1 yn 
对 原始 数据 做 一 次 累加 ,有 
XY = {zDD ,rd (n))}, (2) 


[3 
其 中 zk) 一 Dr Dk = 


由 一 阶 生成 模块 X'" 建立 模型 GM(1,1) ,对 应 的 白化 微分 方程 为 
此 了 名 


十 ar =6, (3) 
其 中 a 和 8 为 待 辨识 常数 . 其 最 小 二 乘 解 为 
$= (a,6) = (BTB)-1 BrY， (4) 
其 中 
Y = [x (2) ,7 (n)]", 
zm(2) 1 
| $s 
—zV(n) 1 


zk 十 D) 二 二 [zDD 十 zk 十 D] 为 背景 值 
方程 (3) 的 离散 解 为 


5 十 D 一 


Rils 
中 
而 
a 
al|e 


(5) 


还 原 的 原始 数据 为 
Fo (十 1D) 一 二 D( 十 1) 一 FPC) 


=(1 一 ez -Ee 
由 公式 (4) 可 知 , 拟 合 和 预测 精度 取决 于 常数 a 和 65, 而 a 和 6 的 求解 则 依 
赖 于 背景 值 z”(& 十 1) ,背景 值 >”(k 十 1) 的 值 成 为 直接 影响 GM(1,1) 模 型 精 
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度 和 适应 性 的 关键 因素 . 


9.2 GM(1,1) 模 型 背景 值 的 改进 


前 面 背景 值 的 求法 实际 上 就 是 数值 积分 中 的 梯形 公式 ,而 梯形 法 的 误差 较 
大 ,精度 较 低 ,为 此 我 们 提出 用 基于 连 分 式 理论 的 有 理 插值 与 广义 梯形 公式 来 
重 构 背景 值 . 

定义 9.125 设 {a,),{b} 为 两 个 实数 列 , 称 形 如 


加 十 一 (6) 
己 十 全 
一 一 Qs 
aa 
所 十 页 十 于 
的 分 式 为 连 分 式 (continued fractions), 记 作 
ba + Kan/b,). 
而 式 
,Pe 3 (7) 
bi 十 2 加 
bz 十 外 
六 十 一 一 全 一 一 
Fe 各 


称 为 连 分 式 (6) 的 nn 次 浙 近 连 分 式 . 其 运算 法 则 按 一 般 连 分 式 运算 , 见 文献 [2]. 
定义 9. 255 下 述 形式 的 连 分 式 : 


bb 十 XT—xo (8) 
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定义 9.30 设 X= {zmvzr…yz，…} 是 实 平面 上 一 点 集 ，F(z) 是 定义 
在 G 字 X 上 的 函数 , 令 
gp[zi] = f(z),i= 0,1,2,.., 


?Im PTR 


TIT 
9 Lz Zio x — plz rT 


gp [zi Tj Ti = 
称 由 上 述 公式 确定 的 9p[zoyzi…zi] 为 函数 7Cz) 在 点 zoyzi，…zi 处 的 1 阶 


逆差 商 . 
定理 9.12 设 


R,(7) 一 9[z] 十 一 一 
9p[m,z] 十 2 
CT A a 


ZX— Zr 
MTT 


(9) 
其 中 p[zoyz，…zi] 天 0,co, 一 0,1，m 为 FCz) 在 zzi…zt 处 的 R 阶 
逆差 商 , 则 有 
及 (zi) = f(zi),i=0,1,,n, 
即 函数 R,(z) 为 函数 f(z) 在 点 zyzi，…yzw 处 的 有 理 插值 函数 . 


将 积分 区 间 [a,6] 划分 为 等 份 , 步 长 hn 二 和 2 二 4 , 则 广义 梯形 公式 为 


m 
Tf) 二 [和 fC@) 十 fa 二 有 D 十 十 f(a 二 (m 一 Dh) 十 去/(O)] 


基于 广义 梯形 公式 的 背景 值 改 进 法 的 步骤 如 下 : 
(1) 取 (2? 式 中 的 一 次 累加 序列 : 
XD = {zDD ,zn)). 
(2) 取 y(k) 二 有 k= 二 1,2,…,n ,Mm 二 4( 或 m= 二 8). 
(3) 把 [y(k) ,x(k)] ,一 1,…,n ,作为 对 应 曲线 上 的 等 间距 点 的 坐标 ， 
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用 有 理 插值 函数 R, (xz) 求 出 以 下 横 坐 标 


3(k 二 二),y(A+ 本) (十 量 ) 大 一 12 一 1 
( 当 m 二 4 时) 对 应 点 的 函数 值 分 别 为 
(kt 二 到) ,zx (8 二 二) ,x(k 二) 二 152 一 1 
(4) 构 造 背景 值 : 
当 m 二 4 时 ,有 


(十 DD = 于 [于 ze 十 zo (8 十 十) 二 zm (k 二 机) 
十 zu (k 十 吝 )+ 生 zk 二 DJ， 
当 m = 二 8 时 ,有 


zz (k 十 1) 


oo 
一 
cl 一 
和 
加 
和 
一 
EF. 
站 


) 

Hzo (hk 二 总) 二 zm (4 二 二)+zm(k 十 宫 ) 
3 7 
8 


zn (kt 主 )+z (t+ 天 诸 寺 zt+D]， 


4 


za (k 二 1) 3 (k+l) $k Fs 


9.3 实例 分 析 


【实例 1】 以 我 国人 均 能 源 消 耗 量 的 预测 为 例 , 比较 本 文 与 文献 [5] 模 型 的 
模拟 预测 精度 ,数据 来 源 于 《中国 统计 年 鉴 》 以 1998 一 2004 年 的 数据 建 模 , 预 
测 2005,2006,2007 年 的 数据 ,结果 如 表 9. 1 所 示 . 按 本 文 方法 ,建立 我 国人 均 
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能 源 消耗 量 灰色 预测 模型 如 (10) 式 所 示 ,文献 [5] 中 所 模型 如 (11) 式 所 示 ，: 
全 9 (十 1) = 1789. 526374eto osrxb 一 1673. 626374， 


和 oO (十 1) = 110. 4380374eo 4 ,k= 1,2,°,n, (10) 
全 (1) = 115. 9; 
全 @ (十 1) = 1792. 315094etooesexe 一 1676. 415094， 
全 o (十 1) = 110. 4419604eto oaxb ,一 1，2， sm (11) 
人 (1) = 115. 9. 
表 9. 1 我 国人 均 能 源 消耗 量 预测 比较 
a 本 文 方法 | 文献 [5] 方 法 
年 份 序号 zol(A) 模拟 数据 ”| 相对 误差 | ”模拟 数据 『 相对 误差 
| | aml | (%) A (k) (%) 
1998 | 1 115.9 115. 90 0 115. 90 0 
| 
1999 | 2 121.4 117.70 3.05 117. 69 | 3.05 
| | 
2000 | 3 126. 4 125. 44 0.76 125. 42 0.77 
200 ee 130.3 133. 69 一 2.60 133. 66 —2.58 
2002 | 5 136. 9 142. 49 一 4.08 局 142. 44 一 4.04 
2003 | 6 153.9 | 151. 86 1. 33 151.79 1.37 
2004 | 7 164. 2 161. 85 人 1.43 161.76 1.49 
2005 
预测 8 179.9 172. 49 4. 12 172. 38 4. 18 
2006 这 
预测 | 9 194.7 183. 84 5.58 183.70 5.65 
| 
2007 
(预测) 10 203. 3 195. 93 3. 63 195. 76 3.71 


将 所 建 模型 进行 平均 相对 误差 和 后 验 差 检验 ,本 文 所 提出 改进 模型 的 平均 
相对 误差 为 1. 46% ,文献 [5] 所 提出 改进 模型 的 平均 相对 误差 为 1. 52% ;计算 
方差 比 二 S/S; ,其 中 S? 为 原始 序列 的 方差 ,S? 为 残 差 序列 的 方差 ,经 计算 ， 
本 文 提出 的 方法 所 建立 的 模型 的 方差 比 为 0. 024 ,而 文献 [5] 提 出 的 模型 的 方差 
比 为 0. 025. 小 误差 概率 P 二 P{|e(k) 一 引 < 0. 6745S; } ,经 计算 ,本 文 的 小 误 
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差 概率 为 1, 所 以 本 文 所 建 模型 的 后 验 差 比 文献 [5] 的 后 验 差 还 要 小 . 综 上 所 述 ， 
应 用 本 文 所 建立 的 模型 具有 更 高 模拟 预测 精度 , 优 于 文献 [5] 的 方法 ,因此 适用 
于 预测 我 国人 均 能 源 消 耗 量 . 

【实例 2】 以 我 国航 空 客运 量 预 测 为 例 ,比较 本 文 与 文献 [5] 的 模型 模拟 预 
测 精度 ,数据 来 源 于 (中国 统计 年 鉴 》 以 2000 一 2004 年 的 数据 建 模 ,预测 2005， 
2006,2007 年 的 数据 ,结果 如 表 9. 2 所 示 . 按 本 文 方法 建立 我 国航 空 客运 量 灰色 
预测 模型 如 (12) 式 所 示 ,文献 [5] 中 所 模型 如 (13) 所 示 : 
人 Wk 十 1) = 41583. 8039e'® sp 一 34861. 8039， 
Vk+1) = 6116. 517487e® sex ,k= 1,2,%,n, (12) 
\ 和 ee(1) = 6722; 
全 (R 十 1) = 42980. 6835ed st5xe 一 36258. 6835， 


全 0 (R 十 1) = 6152. 964555e 1545xb ,k= 1 ,2 (13) 
人 (1) = 6722. 
表 9.2 我 国航 空 客运 量 预测 


本 文 方法 文献 [5] 方 法 
模拟 数据 ”| 相对 误差 | ”模拟 数据 | 相对 误差 
0 (k) (%) A (k) (%) 
6722. 00 0.00 6722. 00 0.00 


7180. 97 —4.56 7171. 34 一 4.69 


8380. 72 8408. 08 一 2. 16 


2003 8759 9780. 92 11.67 . 9858. 09 12. 55 


2004 | 5 12123 11415. 07 一 5.84 11558. 17 一 4.66 
2005 

(预测 ) 6 13827 13322. 23 | 3.65 13551. 44 L090 
2006 

(预测 ) 党 15967 15548. 03 一 2.63 15888. 46 一 0.50 
2007 


8 18576 18145. 71 一 2. 32 18628. 50 0. 28 
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将 所 建 模 型 进行 平均 相对 误差 和 后 验 差 检验 , 本文 所 提出 改进 模型 的 平均 
相对 误差 为 一 0. 17%, 文 献 [5] 所 提出 改进 模型 的 平均 相对 误差 为 一 1. 40%; 计 
算 方 差 比 c = S/S? ,其 中 S? 为 原始 序列 的 方差 ，S; 为 残 差 序列 的 方差 ,经 计 
算 , 本 文 提出 的 方法 所 建立 的 模型 的 方差 比 为 0. 013, 而 文献 [5] 提 出 的 模型 的 
方差 比 为 0.01; 小 误差 概率 P 二 P{|e(k) 一 a| 二 0.6745S?) ,经 计算 ,本 文 所 
用 方法 的 小 误差 概率 为 1, 所 以 本 文 所 建 模 型 的 后 验 差 比 文献 [5] 的 后 验 差 还 要 
小 . 综 上 所 述 , 应 用 本 文 所 建立 模型 具有 比较 好 的 模型 精度 , 优 于 文献 [5] 提 出 
的 方法 ,可 以 用 于 预测 我 国航 空 客运 量 . 
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基于 向 量 连 分 式 理论 的 MGM(1,n) 模 型 


10.1 多 变量 灰色 MGM(1,n) 模 型 


灰色 模型 是 灰色 系统 理论 的 一 个 重要 内 容 , 自 邓 聚 龙 教授 首次 提出 的 20 
多 年 来 ,灰色 模型 在 经 济 管理 等 众多 领域 得 到 了 广泛 的 应 用 "1, 许多 学 者 对 
GM(1,1) 模 型 进行 了 广泛 的 研究 ,为 了 提高 模型 的 拟 合 精度 和 预测 精度 ,提出 
了 一 些 改进 措施 “1. 然而 GM(1,1) 模 型 仅 利用 单一 序列 模拟 和 预测 ,没有 反 
映 多 个 相互 影响 、 相 互 作用 的 变量 之 间 协 调 发 展 和 制约 的 情况 ,而 GM(1,n) 模 
型 反映 了 n 一 1 个 相关 因素 变量 序列 对 系统 特征 序列 的 一 阶 导数 的 影响 ,反映 
了 系统 特征 序列 的 变化 规律 ,如 果 相 关 因 素 变 量 序列 的 预测 值 没有 求 出 ,那么 
系统 特征 序列 的 预测 值 就 不 能 求 出 . 在 实际 的 社会 .经济 系统 中 ,往往 包含 多 个 
变量 ,这 些 变 量 是 相互 关联 、 共 同 发 展 的 ,每 一 变量 的 发 展 变化 都 不 是 孤立 的 ， 
一 个 变量 要 受到 其 他 变量 的 影响 ,同时 也 影响 着 其 他 变量 . 所 以 ,多 变量 灰色 模 
型 (multi-variable grey model) ,简称 MGM(1,n) 模 型 ,从 系统 的 角度 对 各 变量 
进行 统一 描述 . 多 变量 灰色 MGM(1,n) 模 型 是 GM(1,1) 模 型 在 元 变量 情况 
下 的 自然 推广 ,但 不 是 GM(1,1) 模 型 的 简单 组 合 ,也 不 同 于 GM(1,n) 模 型 只 建 
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立 单个 一 阶 微分 方程 ,而 是 建立 个 一 阶 微分 方程 . 

文献 [6] 提 出 了 多 变量 灰色 MGM(1,n) 模 型 ,通过 实例 对 比 说 明了 
MGM(1,n) 模 型 和 GM(1,1) 模 型 ,文献 [7,8] 成 功 应 用 了 MGM(1,n) 模 型 . 然 
而 在 实际 应 用 中 , MGM (1,n) 模 型 存在 一 些 预测 偏差 过 大 的 情况 ,反映 了 
MGM(1,n) 模 型 的 实用 性 有 待 提高 . 本 文 利用 连 分 式 理论 建立 多 变量 灰色 
MGM(1,n) 模 型 的 背景 值 的 计算 公式 ,进一步 提出 了 处 理 多 变量 数据 序列 的 灰 
色 模 型 ,扩展 了 多 变量 灰色 MGM(1,n) 模 型 的 适应 范围 . 最 后 通过 一 个 实例 , 建 
立 中 国 国有 单位 和 城镇 集体 单位 的 多 变量 灰色 模型 . 

在 实际 的 社会 经 济 系统 和 工程 系统 中 ,一 般 都 含有 多 个 变量 , 且 各 变量 
间 是 相互 影响 .相互 关联 ,共同 发 展 的 . 每 一 个 变量 都 受到 其 他 变量 影响 , 同 
时 也 影响 着 其 他 变量 . 多 变量 灰色 MGM(1,n) 模 型 能 够 较 好 地 反映 系统 中 各 
变量 之 间 相互 影响 、 相 互 制约 的 关系 . 当 n = 1 时 ,模型 就 退化 为 GM(1,1) 
模型 . 

令 TI (k) ,i 二 1,2,…sn, 为 第 i 个 灰 时 间 序 列 ; zf (CA) ,i 二 1,2,%… ,nn 
为 相应 的 1 一 AGO, 即 一 次 累加 生成 序列 ， 


4 
XR) = > zfo ,k= 1,2,%%,m. (1) 
j=] 
MGM(1,n) 模 型 的 白化 微分 方程 组 形式 为 : 
0 
二 =a 十 az 入 十 一 十 qz 四 十 权 ， 
on 
aazf 十 azzz 和 2 十 … 十 azvzo 十 ba， 
中 
和 一 az 人 十 az 十 … 十 az 十 加 


为 了 书写 方便 ,写成 矩阵 形式 


YT =AXw +B, (2) 
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其 中 


A= 


dll 3 | 
CR 
dn ”adm 


XY = (zf 0) ,x (0) ,0 x 2))T, 
B= (b ,bb,)T. 
式 (2) 两 边 左 乘 e* , 取 [1, 习 积分 ,得 到 连续 时 间 响 应 为 : 


XP) 一 PKVD+AT eT — 1 B, (3) 
其 中 矩阵 指数 函数 
a SA 
=I1++ 生 二 + 哇 宛 所 


10.2 MGM(1,n) 模 型 背景 值 的 改进 


1. MGM(1,n) 模 型 的 背景 值 
为 了 辨识 参数 A 和 B ,将 式 (2) 离 散 化 ,以 I 个 (十 1) 一 x1? Ck) 为 第 i 个 灰 
导数 向 量 中 的 第 个 分 量 ,而 以 zt? (Ck 十 1) 为 其 对 应 的 背景 值 . 文献 [1] 中 
zf (k 二 1) 一 去 GzipCD 二 zk+D) ,于 是 由 第 i 个 白化 微分 方程 得 到 第 i 个 
灰 微 分 方程 
zk 十 1) 一 Do lt + DD], 


i=1,2,°% ,nk = 1,2,. ,mCO— 1, (4) 
2. MGM(1,n) 模 型 背景 值 的 改进 


定义 1 四 设 {a,) 和 {6b) 为 两 个 实数 列 , 称 形 如 


一 一 一 一 一 (5) 


2 十 一 一 一 
4 
名 十 万 全 


的 分 式 为 连 分 式 , 记 作 如 十 区 (a/6,) ,而 式 


信守 一 = 一 和 一 (6) 


为 连 分 式 的 n 次 渐进 分 式 . 


d 
定义 2 设 v= (wv,…,w) 是 一 个 d 维 实 向 量 , 1v| 二 (Dv?*)' 是 
j=l 
向 量 的 模 , 则 定义 广义 逆 如 下 : 


a 人) 


定义 3 外 设 v 一 (ww，%…uw) 是 一 个 dd 维 实 向 量 , 令 
9[zi] 一 ww 一 0,1… 


= Te 一 Ze 
?Fn 一 FT 二， 


gp [zi ,Tis Th 二 


pp [zi ToT — 9 Lx, zi] 
称 由 上 式 确定 的 [ze ,zi，…zi] 为 关于 向 量 集 V" 在 zo ,xz ，…,z 处 的 ! 阶 向 
量 逆差 商 . 


定理 1 设 
Rlz) = pLm]+ 人 
9[mon] 十 
AER 十 


其 中 p[Lzo,z，…zo] 天 0,co, 一 0,1,…,m, 为 关于 向 量 集 V” 在 ro ,zi ，… ,x 
处 的 & 阶 逆差 商 , 则 有 : 
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R(x) = v= (x Dz DD, 7 0)) ,i = 0,1,° 0n. 


将 积分 区 间 [4,5] 划分 成 m 等 份 , 步 长 h 二 (5 一 a)/m , 则 广义 梯形 公式 为 


T,(/) = 4 十 Fa 十 月 十 … 十 Fa 十 (mm 一 1 有 十 寺 FO) 

当 m 二 4 时 ， 

sj 二 D 二 主 [ 吉 (十 … 十 xz (k 十 总) 二 去 zk+D] 

当 m = 二 8 时 ， 

a = 1fljw rT Re 
ze 十 1D 一 井 [ 羡 zo( 十 + 可 (+D]. 

组 合 公 式 为 : 

(十 D) = 入 eV 二 DD 一 二 sf (二 Dk = 1 ml 

以 2 (十) 一 条 se 中 (Ck 十 1D 一 二 fCk 十 1) 作为 灰 导 数 向 量 的 背景 人 


对 于 本 文 的 向 量 广义 梯形 公式 以 及 外 推 公式 ,可 分 别 对 其 每 一 分 量 进行 即 
可 . 记 
Qi = (daraa da sbi) ,i = 1,2 ,n. 
由 最 小 二 乘法 得 到 a; 的 估计 值 公 为 : 


= (A ,ds db) = LTTILTY,i = 1,2,000n. 


其 中 


竺 oefn (2) raf (2) 1 
3 


六 xzefb(3) + 地 x zffo(3) … 1 


局 
ll 
oo 


年 * em) 十 二 x Dm | 


Yi= (2 (2) ,2 (3) ,zf Cm))T. 
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则 得 A 和 B 的 估计 值 信和 名, 即 


MGM(1,n) 模 型 的 计算 值 为 : 

KH) = eA XPD 4A estan 一 D .Bg = 1,2,°..,m. 
Yo) = XOD), 

OR)= ND NPD) 


一 (IT 一 外 Deheakeh 一 DCXo (1 一 全 一 2,3,.,m. 


10.3 实例 分 析 


为 了 简便 ,这 里 选用 两 个 序列 ,以 1998 一 2005 年 中 国 国有 单位 在 岗 职 工人 
数 作为 序列 zf”, 且 以 1998 一 2005 年 中 国 城镇 集体 单位 在 岗 职 工人 数 为 序列 
zx8) ,建立 zf 和 zi 的 MGM(1,2) 模 型 . 同时 按照 文献 [8] 的 方法 建立 
MGM(1,2) 模 型 ,并 且 分 别 建立 zf? 和 xz 的 GM(1,1) 模 型 .用 以 上 三 种 模型 
模拟 1998 一 2005 年 的 实际 值 ,并 给 出 2006 年 的 预测 值 ,从 而 通过 比较 来 检验 
模型 的 精度 . 

对 序列 xf? 作 光 滑 性 检验 ,得 到 当 & 二 3 时 ,光滑 比 o(&) E (0,0.5) , 且 随 
着 4 的 增 大 , p(k) 递减 . 所 以 zf 为 非 负 准 光滑 序列 . 又 由 文献 [1] 中 的 定理 知 ， 
zf 的 一 次 累加 生成 序列 xf” 具有 准 指数 规律 . 同 理 xz? 的 一 次 累加 生成 序列 
z 刀 也 具有 准 指数 规律 . 
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两 序列 的 MGM(1,2) 模 型 为 : 
一 0.0403zf 一 0. 5741zg2 十 9431. 2969， 
所 = 0.0168zf 一 0. 2388zg2 十 2098. 179. 
文献 [5] 的 MGM(1,2) 模 型 为 : 
oe = 0. 0408z{? 一 0. 5763xg2) 十 9422, 361， 
时 一 0.017zf 一 0. 2394zg2 十 2094. 1965. 
序列 zf? 的 GM(1,1) 模 型 为 : 
3 = 0.0507zf 十 8903. 2867. 
序列 zi 的 GM(1,1) 模 型 为 : 
4 = 0.1328zg2 十 2001. 1119. 


记 文 献 [8] 的 多 变量 灰色 MGM(1,n) 模 型 为 模型 ,本 节 所 建立 的 多 变量 
灰色 MGM(1,nn) 模 型 为 模型 卫 . GM(1.1) ,模型 和 模型 卫 的 模拟 平均 相对 误 
差 和 预测 误差 , 见 表 10. 1 和 表 10. 2. 

表 10.1 三 种 模型 对 zf 的 模拟 值 及 相对 误差 


GM(1,1) 模型 I 模型 下 
序号 | 属 娩 序列 | 模拟 数 据 | 相对 站 天 | 村 拟 数据 | 相对 交差 | 楼 抽 数据 | 相对 训 

1 8809 | 8809 0.000 8809 0. 000 8809 0. 000 
2 8336 8245 1.087 8370 0.408 8378 0. 504 
3 7878 7838 0.514 7813 0. 825 7819 0. 749 
4 7409 | 7450 0. 550 7356 0.715 7361 0.648 
5 6924 7081 2. 271 6982 0.838 6985 0. 881 
6 6621 6731 1. 660 6674 0.800 6677 0. 846 
7 6438 6398 0. 622 6422 0.249 6425 0. 202 
8 6232 6081 2.416 6215 | 0.273 6217 | 0.241 

平均 相对 误差 (%) 1.140 0. 513 0. 509 


第 10 章 基于 向 量 连 分 式 理论 的 MGM ( 1，# ) 模型 大 


表 19.2 三 种 模型 对 xz) 的 模拟 值 及 相对 误差 


GMI(1,1) 模型 I 模型 工 
模拟 数据 IE 模拟 数据 | 相对 洋基 横 夫 开 所 | he 
1900 0.000 1900 0.000 1900 0. 000 
1638 0. 872 1658 0.363 1660 0.484 
1434 0.904 1426 1.451 1429 1.244 
all 
1256 1.174 1238 0.242 1238 0. 242 
5 1071 1099 2.652 1082 1.027 1082 1.027 
6 951 963 1.226 955 0.421 954 0.315 
7 851 843 | 0.948 | 850 0.118 849 0.235 
8 769 738 4.020 765 0.520 764 0.650 
平均 相对 误差 (%) 1.475 0.518 0. 525 


注 : 表 中 的 平均 相对 误差 为 其 百分数 的 绝对 值 . 


从 表 10.1 和 表 10. 2 可 以 看 出 ,本 节 提 出 的 MGM(1,2) 模 型 拟 合 xf" 和 
zg 的 平均 相对 误差 分 别 为 0. 509% 和 0. 525% ,文献 [8] 提 出 的 MGM(1,2) 模 
型 拟 合 zf? 和 zx 名 的 平均 相对 误差 分 别 为 0. 513% 和 0. 518% ,而 本 节 GM(1， 
1) 模 型 对 zf?” 和 xz 名 拟 合 的 平均 相对 误差 分 别 为 1. 14%, 和 1.475%. 可 见 , 本 
节 提出 的 MGM(1,2) 模 型 和 文献 [8] 提 出 的 MGM(1,2) 模 型 对 于 两 变量 都 有 
很 好 的 拟 合 精度 ,都 优 于 GM(1,1) 模 型 . 

用 得 到 的 三 种 模型 预测 2006 年 的 国有 单位 在 岗 职 工人 数 和 城镇 集体 单位 
在 岗 职 工人 数 ,预测 结果 见 表 10. 3 和 表 10. 4. 

表 10.3 三 种 模型 对 zf 的 预测 值 及 相对 误差 


cm) | 模型 模型 1 | 
年 份 | 实际 值 相对 误差 | 相对 误差 相对 误差 
预测 值 (%) 预测 值 (%) 预测 值 (%) 
2006 | 6170 5781 | 6.312 | 6045 | 2.026 中 6046 | 2.010 
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表 10.4 三 种 模型 对 zj) 的 预测 值 及 相对 误差 


| enD) 模型 I 模型 了 
年 份 | 实际 值 相对 误差 相对 误差 相对 误差 
| 预测 值 | 太守 二 郑 | 预测 值 | 和 全 郑 | 预测 值 | 相合 
2006 | 726 | 646 10.98 | 694 4.408 692 4.683 | 


从 表 10.3 和 10.4 可 以 看 出 ,本 节 提 出 的 MGM(1,2) 模 型 预测 zj" 和 xz 名 
的 相对 误差 为 分 别 2. 01% 和 4. 683%, 文 献 [8] 提 出 的 MGM(1,2) 模 型 预测 
zf2 和 xz" 的 相对 误差 分 别 为 2. 026% 和 4. 408% ,而 本 节 GM(1,1) 模 型 对 zf 
和 xz 预测 的 相对 误差 分 别 为 6. 312% 和 10. 98%. 本 节 提 出 的 MGM(1,2) 模 
型 和 文献 [8] 提 出 的 MGM(1,2) 模 型 预测 精度 都 比较 高 . 因此 ,本 节 提 出 的 
MGM(1,2) 模 型 比较 好 地 反映 了 国有 单位 在 岗 职 工人 数 和 城镇 集体 单位 在 岗 
职工 人 数 之 间 相 互 制约 、 共 同 发 展 的 关系 . 
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非 等 间距 GM(1,1) 模 型 时 间 响 应 
函数 的 优化 


11.1 非 等 间距 GM(1,1) 模 型 的 建 模 机 理 


灰色 系统 理论 自 邓 聚 龙 教授 创立 以 来 ,在 许多 领域 得 到 了 广泛 应 用 ,基于 
贫 信 息 的 灰色 预测 成 功 地 解决 了 许多 信息 不 完全 的 预测 问题 中 但 是 ,灰色 系 
统 模 型 大 都 建立 在 等 间距 序列 的 基础 上 ,而 在 现实 中 得 到 的 原始 数据 可 能 是 非 
等 间距 序列 . 所 以 ,建立 非 等 间距 序列 的 GM(1,1) 模 型 具有 现实 意义 . 许多 学 
者 对 此 进行 了 深入 研究 ,文献 L[10 一 16] 阑 述 了 几 种 非 等 间距 序列 GM(1,1) 模 
型 的 建 模 机 理 ,并 在 具体 的 领域 得 到 了 成 功 应 用 . 但 是 所 建 模型 的 精度 较 低 , 而 
且 还 没有 完全 脱离 等 间距 序列 的 建 模 思想 ,应 用 范围 受到 一 定 的 限制 . 

基于 上 述 局 限 性 ,本 章 将 通过 研究 非 等 间距 原始 序列 的 特点 和 GM(1,1) 
模型 的 建 模 机 理 , 以 原始 序列 一 次 累加 序列 的 观测 值 与 模拟 值 的 残 差 平方 和 最 
小 为 条 件 ,构建 非 等 间距 GM(1,1) 模 型 的 时 间 响 应 函数 的 优化 模型 ,并 将 该 模 
型 应 用 于 文献 [11] 中 钛 合金 疲劳 强度 随 温度 变化 关系 的 研究 ,结果 表明 ,该 模 
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型 在 模拟 精度 和 预测 精度 方面 都 取得 了 较 好 的 效果 . 
定义 11.19" 设 原始 序列 Xo 二 (zx (1) ,x(k2),…,z0(k,)) , 相 邻 


分 量 之 间 的 间 工 为 

Mk; = ki;— ki = 2,3,.°°,n, (1) 
车 从; 为 常数 , 则 称 序列 X 为 等 间距 序列 , 若 从 不 为 常数 , 则 称 Xe 为 非 等 
间距 序列 . 


在 实际 应 用 中 存在 大 量 的 非 等 间距 序列 的 拟 合 和 预测 问题 ,特别 是 在 工程 
技术 领域 , Ak; 不 为 常数 . 为 了 综合 考虑 序列 的 间距 从 ,在 对 原始 序列 进行 一 
次 累加 时 ,以 序列 的 间距 作为 乘 子 ,建立 非 等 间距 GM(1,1) 模 型 . 

定义 11.2 设 序列 Xo = (zo CD)，zo()，…zo (Cn)) ,车 其 中 每 个 
分 量 都 满足 

xD (kh) 一 Deo dasi = 1,2,° on, (2) 
规定 x 中 (ki) 一 x(k) , 则 称 XX 为 XW 的 一 次 累加 生成 序列 . 
定理 11.10" 设 X% 为 非 负 光滑 序列 , 即 
Xe = (xO ki) sz Ck) ,eT ))， 
则 称 
XY = (Vk) ,ID ke) se, TY )) 
为 XW”: 的 一 次 累加 生成 序列 时 ,对 X"” 建立 的 白化 微分 方程 


Etartd() =6 (3) 


的 灰色 微分 方程 为 ; 

TO (Rh D)ARr + az Ck) — BAkin ， (4) 
其 中 zx CH) 为 z(z) 在 [ki,km] 上 的 背景 值 . 灰色 微分 方程 的 最 小 二 乘 参 
数 估计 为 : 


4 = (4,0)7 = (BB)-B'Y, (5) 
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其 中 
TO (ka)Akz —zV ks) As 
TO (ks)Aks —zV (ks) As 
TO (k,)Ak, 一 zV(k) 名， 


11.2 非 等 间距 GM(1, 1) 模 型 时 间 响 应 函数 的 优化 


文献 [12,13] 以 序列 中 第 一 个 分 量 作为 非 等 间距 GM(1,1) 模 型 的 初始 条 
件 ,没有 对 白化 微分 方程 的 时 间 响应 函数 进行 优化 ,不 能 保证 模拟 序列 和 原始 
序列 的 最 优 拟 合 , 造 成 模拟 精度 和 预测 精度 都 比较 低 . 为 了 改进 以 上 问题 ,笔者 
运用 最 小 二 乘法 来 优化 非 等 间距 GM(1,1) 模 型 的 时 间 响 应 函数 . 
定理 11.2 设 Y,B ,人 如 定理 11.1 所 述 , 人 一 (9)7 一 (BTB)-:Bry， 
则 
(1) 和 白化 微分 方程 了 二 2 十 ar (2) 一 的 时 间 响应 函数 为 
六 (ze 之 )e= 
ZD(i) 一 妈 ! 5 e< 十 
De 
:=1 
(2) 非 等 间距 GM(1 ,1) 模 型 zx (RH )A 十 az (CA ) 一 AAA 的 时 间 响 
应 序列 为 


了 


a 


(6) 


D(x) 2) 
z0 kn) 一 也 em 十 也 (7) 


-- a 
3 2 
1 


二 灰色 缓冲 算 子 理论 及 其 用 


(3) 还 原 值 
(kn)= 和 Ca 一 人 2 
1 
zk) 一 也 -a 
四 Er >( ): 和 (8) 


ez 


证 明 :1) 白 化 微分 方程 os 十 ar (1) 一 4 的 通 解 为 ， 


b 


IV) 一 ce 让 
为 了 求 出 c ,构造 一 次 累加 生成 序列 的 残 差 平方 和 函数 RSS(c) : 


RSS(c)= D(z) — "0)) 
‘St 


(i ee dN 
= 2 (zx 4 ) 。 
这 是 一 个 二 次 函数 ,一 定 存在 极 小 值 . 令 区 (人 = 0 , 求 出 唯一 驻 点 ,得 


bp zm —L)e 
全 2 0) (9) 
De 
根据 实际 情况 , 知 函数 RSS(c) 在 此 驻 点 取得 最 小 值 . 
(2) 令 1 二 ,由 (9) 式 可 得 : 
Dr?) 2) 
XIV kin) = SR , 咎 才 
De 人 
(3) 由 累 减 还 原 可 得 . 
定理 11.3 设 原 始 序列 Xo 一 (zt (RD ,zeo Ce) ze (As)), 非 等 间距 
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GM(1,1) 模 型 的 时 间 响 应 序列 累 减 还 原 值 为 


入 GO) A 
wk) 一 芋 Ce Ce) ， 
有 

(1) 当 态 委 于 时 ， 全 o (k;) 为 模型 的 模拟 值 ; 


(2) 当 之 nn 时 , 4 (k,) 为 模型 的 预测 值 . 


11.3 实例 分 析 


为 了 便于 比较 ,采用 文献 [11] 中 的 工程 技术 实例 ,该 序列 是 P. G 福 雷 斯 研 
究 钛 合金 疲劳 强度 随 温度 变化 的 实验 结果 ,这 是 一 个 非 等 间距 序列 . 笔者 对 原 
始 数据 序列 按照 本 文 的 方法 和 文献 [16] 的 方法 同时 建 模 , 然 后 对 比分 析 . 
本 文 所 建立 的 非 等 间距 GM(1,1) 模 型 简 记 为 优化 模型 ; 
AV (kin) 一 一 6. 32371 X 10e ®t 十 5.73646 X 105， 
其 中 i = 1,2,…,8. 
文献 [11] 的 模型 为 : 
全 0 (Ar ) 一 一 5. 77025 X 105e oct 1o0) 十 5. 77585 X 105, 
其 中 i 二 1,2,…,8. 
文献 [16] 的 模型 为 : 
AV (kn) 一 一 5. 8614 X 105e m602sckn 100) 十 5. 867 X 105 ， 
其 中 = 1,2,…,8. 
三 种 模型 的 模拟 值 及 其 相对 误差 见 表 11. 1. 


次 色 缓冲 算 子 理论 及 其 应 用 


表 11.1 三 种 模型 的 拟 合 值 及 其 相对 误差 


文献 [11] 的 模型 文献 [16] 的 模型 优化 模型 
T/C| c-， 序 司 相对 误差 相对 误差 相对 误差 
模拟 值 | “(%) ”| 模拟 值 | 。(%) 。 | 模拟 值 | 。 (96) 
100 560 kl 560 0.00 560. 00 0.00 560. 00 0. 00 
130 | 557.54 | k | 563. 04 0.99 554. 81 0.49 557. 35 0.03 
170 | 536.1 | k | 537.06 0.18 536. 49 0.07 537. 36 0. 24 
210 | 516.1 | A | 512.31 0.73 516. 27 0.03 516. 60 0.10 
240 | 505.6 | k; | 494. 28 2. 24 499. 20 LL27 499. 08 1. 29 
270 | 480.1 | &k; | 481.20 1.01 485. 02 1.03 Et 484. 55 0.93 
310 | 467.4 | k; | 468. 26 0.18 469. 00 | 0. 34 468. 14 0. 16 
340 | 453.8 | ks | 457. 30 0.77 453. 49 0.07 452. 27 0. 34 
平均 相对 误差 (%) 0.76 0.41 0. 39 
注 : 表 中 平均 相对 误差 为 其 绝对 值 
从 表 11. 1 可 以 看 出 ,本 文 提出 的 优化 模型 模拟 精度 (模拟 准确 度 ) 为 
0. 39%, 文 献 [11] 和 [16] 的 模拟 精度 分 别 为 0.76% 和 0.41% ,显然 文献 


[16] 的 模拟 精度 优 于 文献 [11] 的 模拟 精度 ,本 文 所 建 模 型 的 模拟 精度 优 于 文 

献 [16]. 
三 种 模型 的 预测 值 及 其 相对 误差 见 表 11. 2. 

表 11.2 三 种 模型 的 预测 值 及 其 相对 误差 


文献 [11] 的 模型 文献 [16] 的 模型 优化 模型 


| 


相对 误差 相对 误差 相对 误差 
模拟 值 (%) | aa (%) ET (%) 


380 | 436.4 | 


448. 17 一 2.70 | 438. 51 一 0.48 | 436. 95 一 0. 13 


由 表 11. 2 可 知 ,本 文 提 出 的 优化 模型 的 预测 精度 也 高 于 文献 [11,16]. 可 
见 , 文 中 的 优化 模型 较 好 地 反映 了 铁合金 疲劳 强度 和 温度 的 变化 关系 . 
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GMI(1,1I) 模 型 的 病态 问题 研究 


12.1 矩阵 条 件数 


灰色 系统 是 根据 关联 度 、 灰 数 、 灰 导数 等 一 系列 数学 方法 建立 起 来 的 连续 
性 的 微分 方程 , 它 能 利用 少量 的 数据 进行 建 模 . 目前 ,被 广泛 应 用 于 工业 、 农 业 、 
经 济 等 领域 . 在 实际 应 用 过 程 中 ,数据 矩阵 中 的 元 素 一 般 是 观测 得 到 的 ,所 以 不 
可 避免 地 带 有 误差 . 这 对 于 模型 的 参数 辨识 有 一 定 的 影响 . 

定理 12.10 设 AE C0”™,6A € 0”, EC', 车 对 C” 上 的 某 一 矩阵 范 


数 上 .上 有 AT 上 .中 84 1 二 1, 则 非 齐 次 线性 方程 组 Ar 一 0 与 (4 十 
64A)(z 十 民 ) ==b 十 是 的 解 满足 
| sz |。 < 一 4A1x14-1 al ,lol ). 


rl STHAN x IA x (CTAT + TT 


其 中 外， ,是 C" 上 与 矩阵 范 数 ‖。 || 相 容 的 范 数 . 

由 以 上 不 等 式 可 以 看 出 ,数据 的 误差 对 线性 方程 组 解 的 影响 与 Al X 
上 4 中 的 大 小 有 关 . 因此 ,可 以 A 1 X 1 4 一 中 作为 影响 线性 方程 组 解 的 大 
小 的 一 种 度量 . 


天 灰色 缓冲 算 子 理论 及 其 必用 


定义 12.10 设 AEC“% ,| :| 是 C 上 的 矩阵 范 数 , 称 K(A) = 
1A1 关外 4A” | 为 矩阵 A 的 条 件数 . 


当 A 为 正规 矩阵 时 ,有 天 (4) 一 入 ,其 中 ,4, 分 别 是 A 的 模 最 大 和 最 


小 特征 值 . 由 于 A € C*” 的 所 有 甜 阵 范 数 都 是 等 价 的 ,为 了 方便 ,我 们 取 2 
范 数 . 

实践 中 ,一 般 认 为 :车 1 一 K(CA) 二 10, 则 A 为 良 态 的 ;车 10 志 K(A) 二 100， 
则 A 为 轻 度 病态 的 ;车 100 < KCA) 二 1000, 则 A 为 较 强 病态 的 ; 若 1000 < K(A)， 
则 A 为 严重 病态 的 ( 见 文献 [1]). 


12.2 ”GM(1,1) 模 型 的 病态 分 析 


定义 12. 2 设 系统 行为 数据 为 非 负 光滑 序列 Xo 一 (x (1)，…， 
ZV(n)) ,XW 为 XW 的 1 一 AGO 序列 , ZW 为 X'” 紧邻 均值 生成 序列 , 称 
ZW (R) 十 az(k) 二 5 为 灰色 微分 方程 ,也 被 称 为 GM(1,1) 模 型 的 定义 式 . 

定理 12.20] 车 a 二 (a,6)" 为 GM(1,1) 模 型 参数 , 且 


OK(2) 六 [zc (&)]2 和 so 
Y= ,BTIB 一 | 人 ， 
To (n) — Dz (k) n—1 
k=2 


则 灰色 微分 方程 x" (4) 十 az"”(k) = 6 的 最 小 二 乘 估 计 参 数 序列 为 
a= (a,b)T = (BTB)7BY, 


其 中 系数 矩阵 
DY 一 六 xm 
BTB a k=2 k=2 
一 六 oo n—1 
经 


0 
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其 伴随 矩阵 为 


n—1 Dd) 
(BrB)， = Ei 
BDz?) > Lz) 
在 全 
由 定义 12. 1 知 :系数 矩阵 B7B 的 2 范 数 的 条 件数 为 
K(B'B)= || B7ZB 1 x | (CBB)1。 


BB|:. x | (BB)? 2 _ {alx Iar | 
1B87B1 [BB| ， 


其 中 , 1 为 B7B 模 最 大 的 特征 值 , ”为 (B7B)” 模 最 大 的 特征 值 . 
定理 12.30] 设 AE C”, 则 A 的 任 一 特征 值 满足 


[Al<nx max|ai,l. 
3 


由 定理 12.3 得 ， 


TT [lx [27 | 4max|ay |X max|a; | 
KB'B) = a < 1 5 


显然 ,由 B"B 与 (B7B) ”的 表达 式 知 ,矩阵 B7B 与 (B7B) ”的 模 最 大 值 无 
非 就 是 


Dw, nl | Dn)| 
k=2 k=2 
三 者 之 一 . 
引 理 车 a = (41,…sa,) 和 65 二 (b1，,…,b,) 是 实数 序列 , 则 
Tap < (Tat Tn. 
=1 k=1 k=1 
等 号 成 立 的 充 要 条 件 是 序列 a 与 5 线性 相关 . 
由 引 理 知 
(| p37 SD LP Xa. 
下 面 我 们 分 三 种 情形 来 分 析 和 矩阵 B7B 的 条 件数 . 


) 办 色 缓冲 算 子 理 论 及 其 应 用 = 


(1) 对 于 所 有 的 &, 有 zm(k) 二 1,n 一 1 为 矩阵 BTB 与 (BTB)* 中 的 最 大 值 
元 素 , 则 


KCB7B)<< tml er lee lei 
4 (ma 一 1)2 
(一 1) > [zo(b) 卫 一 ( [ze)]) 
k=2 k=2 


4(n—1) 


( 3 [ze C8)])’ 
之 [zo(&) 了 一 一 全 3 


(Dzw 0])’ 
当 3 [zDD 一 一 一 一 之 1 时 ， 
K(B'B) < 4(n—1). 
由 于 灰色 系统 讨论 的 是 小 样本 问题 , x 过 10, 则 K(B7B) 志 40 二 100. 此 时 
显然 GM(1,1) 模 型 不 存在 大 的 病态 问题 . 


(Pzw I) 
当 0 二 2 [zc (A) 了 :一 二 去 1 时 ,条 件数 的 大 小 取决 于 


( 袜 [emco]) 
> [ze (k)]? 一 一 4=2 


= 
接近 于 0 的 程度 . 
(六 [zkD]) 
当 2 [z (了 一 一 一 一 0 时 , 则 
k=2 


(DEev 8])’ = (DF Lv WF) nD. 
=2 k=2 
此 时 ，{z (k)) 近似 为 常数 序列 , 则 {zm (Ck)) 也 近似 为 常数 序列 , 故 {x0 Ck)) 
近似 为 首 项 不 为 0, 其 他 各 项 全 为 0 的 常数 序列 ,对 于 此 序列 建 模 没 有 实际 
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(2) 对 于 所 有 的 上 有 zx (Re) 二 1， > [zn (4)]: 为 矩阵 BTB 与 (BITB) ”中 
的 最 大 值 元 素 , 则 肖 
K(B'™B)< le Bs -| 
4 也 [zo (了 
(一 DD [zo (6) 卫 一 (BE) 


45) Lz? C8) 
二 二 k=2 


(BE? 1) 


(01 
D3 [27 Ck) 
大 一 2 
当 
(> [ze)]) CD LzY Ce)]) 
0 到 -和 委 2 一 2, (一 1 一- 过 1 
> [zc C4) 2 [ze C8) 
k=2 k=2 
时 ， 


K(BIB) < 3 [ze (WF SA4n—1) Jax[z (AD 于 . 
尽管 灰色 系统 研究 的 是 小 样本 问题 ,比如 n 二 5, 但 只 要 <”(A) 一 10, 则 
K(B'B) < a [z® (kJ <4X4X10 = 1600, 


而 矩阵 的 条 件数 K(B7B) 宇 1000, 则 和 矩阵 B7B 就 是 最 严重 的 病态 矩阵 . 实际 问 
题 中 2 (&) 远 远 大 于 10,n 可 以 取 10, 则 病态 问题 肯定 相当 严重 ,文献 [1] 中 的 
算 例 也 表明 这 一 点 . 


A 灰色 缓冲 算 子 理论 及 其 几 用 


CFV 
同样 , 当 n 一 2 过 -与 和 一 一 一 之 n 一 1 时 , 条 件数 同时 取决 于 4》) [zk) 了 了 
> [2 (8) k=2 


CF zm) DD’ 

的 大 小 与 天 -一 一 一 一 接近 于 n 一 1 的 程度 ,比如 当 
by [zu Ck 
k=2 


(DE) 
3 [zm 


—(n—1)= 0.8, 


CP ev) 

此 时 , -全 一 一 接近 于 n 一 1 的 程度 还 是 相当 远 , 但 
> [zcb (kt)]: 
k=2 


K(BT™B) <52 [z™ C8), 
同样 由 于 原始 数据 较 大 ,1 次 累积 和 数据 更 大 ， 比如 zm Ck) 二 100,n 一 5, 从 而 
K(BT'™B) < 5 [zx (k)Y = 5X4X100 = 2x10 > 1000, 
文献 [中 中 的 算 例 显示 逢 阵 B7B 的 的 条 件数 高 达 107 . 
当然 , 若 Dp [ze CA) 了 <50 或 1 思 [z (8) 了 过 100 时, 则 引起 病态 问 


(HEP 
题 的 原因 就 是 -与 一 一 一 一 - 接近 于 n 一 1 的 程度 . 当 
2 Lz" (RF 
CDE AD)’ 
k=2 
2 [z™ (gk) F 
k=2 
此 时 ,{z"(k)}》 近 似 为 常数 序列 , 即 {zm (CA)} 也 近似 为 常数 序列 , 故 {x 中 Ck)} 


一 (2 一 1)， 


@®— 
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近似 为 首 项 不 为 0, 其 他 各 项 全 为 0 的 常数 序列 ,对 于 此 序列 建 模 没 有 实际 


(3) 当 | 六 [zw(b]| 最 大 , 即 
大 一 2 
| 立 [eokp]|> 六 [ob 了 ,| 六 [xbo]|>w 一 1>1， 
k=2 k=2 大 一 2 
则 
CD Een 08) = [Se cI|x [tw]|> 习 [zk)F) nO— 1), 
fa2 t=? t=2 t=2 
由 引 理 知 不 可 能 . 即 | 六 [ze (4)]| 不 可 能 最 大 . 
k=2 
总 之 ,三 个 数 (加 [0]) ,| 六 [ze CD]| ,mn 一 1 中 的 最 大 值 只 能 是 
> 名 
《加 C000》 ,一 1 则 针对 序列 (sw(bD) 中 ,有 的 元 素 大 于 1, 有 的 元 素 小 
于 1, 也 只 能 是 (》 [zw Ck)])* 最 大 或 wn 一 1 最 大 ,可 分 别 归于 (1)、(2) 进 行 
k=2 
讨论 . 
综 上 所 述 , 引 起 矩阵 条 件数 变 大 的 原因 只 有 两 个 ,原始 数据 较 大 或 病态 
数据 ( 即 近似 为 首 项 不 为 0, 其 他 项 为 0 的 常数 序列 ). 对 于 后 一 类 问题 ,没有 
现实 意义 ,不 必 深 究 . 对 于 原始 数据 较 大 ,可 选择 数 乘 变 换 , 比如 对 非 负 原始 


数据 
Xe = (x0(1) ,zn)), 


及 其 一 次 累积 数据 

KY = (zDD ,ez Cn))， 
紧邻 均值 生成 序列 

ZY = (2 (2) ,2 Cn)), 
显然 有 


XPD) SEIV Nn) ,2 (2) < Sz (n). 


) 灰色 缓冲 算 子 理论 及 其 K 用 


对 原始 数据 zx (&) 乘 以 一 个 数 p < -ml ， 则 得 新 的 数据 


Yk) 一 0XzD() LpXz Yn) 1 k= 2 n. 
故 可 转 人 本 文中 的 (1) 来 解决 . 从 而 对 于 新 的 数据 y”(k) 建 模 应 该 不 会 有 太 大 
的 病态 问题 . 对 得 到 的 结果 进行 反 变换 ,文献 [1] 已 详细 讨论 . 

若 1 二 K(A) 二 10, 则 A 为 良 态 的 .车 10 志 K(A) 一 100, 则 4 为 轻 度 病 
态 的 . 车 100 过 K(A) 二 1000, 则 A 为 较 强 病态 的 . 若 1000 之 K(A), 则 A 为 严 
重病 态 的 . 据 此 论述 ,最 后 给 出 一 个 如 何 解决 灰色 GM(1,1) 建 模 中 病态 问题 的 
具体 方案 : 

(1) 若 ?一 1 为 矩阵 B 了 3 与 (B"B)" 中 的 最 大 值 元 素 , 可 直接 对 原始 数据 建 模 . 

(2) 若 > [Lz (&) 耻 为 矩阵 B7B 与 (B87B)” 中 的 最 大 值 元 素 , 则 

D3 [ze C8) 


(I) 当 0 过 到 一 一 
了 [ze Ch) 
k=2 


入 "一 2 时 ， 


车 max[Lz" (k)] < 时 , 1 二 K(BB* ) 一 10, 则 矩阵 BB" 为 良 
态 的 . 可 放心 利用 GM(1,1) 建 模 ; 


车 max[z"(h)] 和 zo; 时 , KCBB*) < 100, 则 矩阵 BB 为 良 态 
的 或 轻 度 病态 的 ,可 以 利用 GM(1,1) 建 模 ; 


若 max[z" (和 )] < Too 时 , KCBB*) 一 1000, 则 矩阵 BB 为 较 强 


病态 的 ,可 以 谨慎 利用 GM(1,1) 建 模 ,最 好 对 原始 数据 先 施行 数 乘 变换 ,然后 
建 模 ; 
车 max[z"(k)] 之 zt; 时 ， K(BB*) 有 可 能 大 于 1000, 则 和 矩阵 


BB" 为 严重 病态 的 ,建议 对 原始 数据 先 施行 数 乘 变换 ,然后 建 模 ; 


@@= 


第 12 章 “GMI(1,1) 模 型 的 病态 问题 研究 多 
(> [ze (8) 


2 二 人 委 ? 一 1 时 ,比如 当 
3D} [x (k)J 
k=2 


CD Ez RT 
0 入 一 1 一 所 :一 一 一 = 入 1， 
> [zx (k)J 


t=2 


且 a 不 太 小 时 ,比如 a 宇 0.5 时 ,可 分 别 根据 
Bex" (EY Bo’ 


Fy) 
BL (的 ] 科 VER 一 TD， 


(CI) 当 ?= 


(CD 1000 
i 04] 4a(Cm 一 1) 


CD 1000 
BL" JV 
采用 类 似 于 ( 工 ) 中 的 方法 进行 建 模 . 
数 乘 变换 中 的 p 选 择 , 只 要 p< zk 即 可 . 


全 
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